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ABSTRACT 

Aims. The distance of NGC 1316, the brightest galaxy in the Fornax cluster, provides an interesting test for the cosmo- 
logical distance scale. First, because Fornax is the second largest cluster of galaxies within S; 25 Mpc after Virgo and, 
in contrast to Virgo, has a small line-of-sight depth; and second, because NGC 1316 is the single galaxy with the largest 
number of detected Type la supernovae (SNe la), giving the opportunity to test the consistency of SNe la distances 
both internally and against other distance indicators. 

Methods. We measure surface brightness fluctuations (SBF) in NGC 1316 from ground- and space-based imaging data. 
The sample provides a homogeneous set of measurements over a wide wavelength interval. The SBF magnitudes, cou- 
pled with empirical and theoretical absolute SBF calibrations, are used to estimate the distance to the galaxy. We also 
present the first B-band SBF measurements of NGC 1316 and use them together with the optical and near-IR SBF 
data to analyze the properties of field stars in the galaxy. 

Results. We obtain (m-M) = 31.59 ± 0.05(stat.) ± 0.14(sys.) mag, or d = 20.8 ± 0.5(stat.) ± 1.5(sys.) Mpc. When 
placed in a consistent Cepheid distance scale, our result agrees with the distances from other indicators. On the other 
hand, our distance is 17% larger than the most recent estimate based on SNe la. Possible explanations for this 
disagreement are the uncertain level of internal extinction, and/or calibration issues. Concerning the stellar population 
analysis, we confirm the results from other spectro-photometric indicators: the field stars in NGC 1316 are dominated 
by a component with roughly solar metallicity and intermediate age. A non-negligible mismatch exists between B-band 
SBF models and data. We confirm that such behavior can be accounted for by an enhanced percentage of hot horizontal 
branch stars. 

Conclusions. Our study of the SBF distance to NGC 1316, and the comparison with distances from other indicators, 
raises some concern about the homogeneity between the calibrations of different indicators. If not properly placed in 
the same reference scale, significant differences can occur, with dramatic impact on the cosmological distance ladder. 
Our results on the stellar populations properties show that SBF data over a broad wavelength interval are an efficient 
means of studying the properties of unresolved systems in peculiar cases like NGC 1316. 

Key words, galaxies: elliptical and lenticular, cD - galaxies: distances and redshift - galaxies: clusters: individual: 
NGC 1316 - galaxies: photometry - galaxies: stellar content - galaxies: peculiar 

1. Introduction range of applicability in terms of distances and of useful 

targets in order to minimize the error propagation over the 
The past three decades have seen remarkable progress cosmological distance scale; and (iii) improve/analyze the 
the study of the distance scale o f the Universe matchmg between mdependent mdicators. A particularly 



(iJacobv et all [1992: Fcrrarese et a l. 2001 iFreedman et all promising distance indicator for addressing the three listed 



I200 H iFreedman fc MadordlMoD . resulting in a gen- hems is the surface brightness fluctuations technique (SBF 

^ Tonrv fc Schneider 1988; Tonrv et al. 1990; Blakesle e etall 

20090 . 



eral convergence of distances based on different indica- 
tors. Nevertheless, to resolve the lingering disc repancies 



(jTammann et al.ll2008l : IPVeedman fc Madore 120101 ) . there is After the first two "rungs" of the cosmological distance 



an urgent need to (?) lower the statistical (intrinsic) and sys- scale, represented by geometric methods and by primary 
tematic (external) errors for single distance indicators; (ii) indicators (variable stars, main-sequence fitting, etc.), the 
obtain distance measurements from indicators with a large SBF method is one of the most accurate indicators, with a 
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median 0.14 mag accuracy on distance moduli, or ~ 7 % in 
dista nce, up to ~ 100 Mpc (|Mei et al.ll2003t iBiscardi etaP 
l2008f) . However, the accuracy gets considerably smaller, to 
a mean ^ 0.08 mag (^^4% in distance) from space-based 
optical data (Table [J) 

By definition, the SBF signal corresponds to the ratio 
of the second to the first moments of the luminosity func- 
tion of stars in a galaxy. As opposed to the surface bright- 
ness that does not scale with the distance of the galaxy, 
the SBF signal scales inversely with the distance squared. 
Observationally, the SBF method relies on the measure- 
ment of the intrinsic flux variance in a galaxy, generated by 
the Poissonian fluctuations in the surface brightness due 
to the statistical variation of the stellar counts in adjacent 
resolution elements. The variance, normalized to the local 
mean surface brightness, is converted to an apparent mag- 
nitude, m, from which the distance modulus, fh—M, follows 
once the absolute M is known. 

Given its definition, M in a given bandpass is dependent 
on the properties of the underlying stellar populations. The 
analysis of large samples of early- type systems, including 
galaxy in groups, has made it possible to characterize the 
dependence of M on stellar population properties using lin- 
ear relations with respect to some b r oad-band optical colou r 
(|Tonrv et al.ll200lt iMei et al.ll2007t ICantiello et "alll2007a[l . 

In Table [1] we report the median errors on SBF mea- 
surements, S{fh), and on the associated distance moduli, 
6 (fh—M), derived from different samples. For optical SBF, 
the typical accuracy for ground-based measurements has 
been ^ 0.2 mag (10% in distance). Because this is dom- 
inated by measurement errors, the superior resolution of 
HST leads to an improved accuracy of ^0.08 mag. For a 
comparison, the internal scatter of the period-luminosity 
(PL) relation of Cepheids ranges from 0.20 mag in the V 
band, to 0.0 9 mag in the Spitz er [3.6] /im and [4.5] fim 
bands I'e.g.. iNgeow et al.ll2009[ ). The mean uncertainties 
for the near-IR SBF sample quoted in Table [T] are gener- 
ally larger than the optical ones, ^ 0.2 mag (~10% on 
distances) even in the case of space-based data. This is 
due both to observational/technical iss ues (e.g. dark cur- 
rent p atterns, "wormy" background, see lJensen et al.l[l998l 
1200 It) , and to the scatter of the calibration (mostly re- 
lated to the sensitivity of the SBF signal to the properties 



of AGB and TP-AGB stars in these bands, see IMei et al 



2001 




2005 





pr 



Liu et al] 20021: Jensen et al] 20031 : iRaimondo et al . 
Gonzalez-Lopezlira et al.ll2010t) . The availability of 



Wide Field Camera 3 (WFC3, on board of the Hubble 
Space Telescope, HST) is expected to significantly improve 
the situation in this wavelength regime. Finally, the zero 
point of the calibration is typically tied t o the Cepheid dis- 
tance scale to an accu racy ^ 0.08 mag (jTonrv et al.ll2000t 
iBlakeslee et"alll2010D . 

Thus, SBF distances are characterized by reasonably 
well defined and small (especially in optical bands) inter- 
nal and calibration errors (item (i) in the list above). The 
technique has been used to estimate distances for Local 
Group galaxies (even closer than that if one t akes into ac- 
count the work on Galactic globular clu sters by Aihar et al.l 
19941) . out to galaxi es at > 100 Mpc (| Jensen et al.ll200lt 



Biscardi et al.ll2"008f) . With the highly improved near-IR 



imaging capabilities of the WFC3/IR and similar instru- 
ments, and thanks to the much brighter SBF signal in the 
near-IR, the upper limit on SBF distances is expected to 
increase significantly. Thus, SBF can encompass more than 



two orders of magnitude in distance, bridging local to cos- 
mological distances with the use of a single indicator (item 
(ii) above). 

In this paper we present measurements of SBF 
magnitudes for the intermediate-age merger remnant 
NG C 1316, al so known as Fornax A (e.g., ISchweizei1ll980l : 
Tc rlevich &: Fo rbes 2002i) . This galaxy is peculiar in many 
ways. Although by far the brightest member of the Fornax 
cluster, it is not near the cluster center, being ~ 3.7° 
away from the central giant elliptical NGC 1399 (projected 
separation of ~ 1.3 Mpc at the distance of Fornax). It 
shows numerous dust features, a prominent dust lane. 
Ha filaments, lo ops, and tidal tails originally analyzed by 
iSchweizeri (fT98l . Moreover, it is a powerful radio-galaxy 
and, to date NGC 1316 is the single galaxy with the largest 
number of discovered Type la supernovae (SNe la hereafter; 
four events recorded). The latter property makes NGC 1316 
a remarkable place to test the extragalactic distance scale, 
because of the role of SNe la for cosmological distances, 
and because NGC 1316 is one of the nearest massive post- 
merger galaxies. 

We have collected data covering a large wavelength in- 
terval (from B to H band), with the specific purpose of 
performing a self-consistent analysis of SBF data for this 
galaxy in order to carry out a comprehensive study of the 
SBF in NGC 1316, and consequently, of the galaxy distance 
(item (iii) above). 

Furthermore, given their dependence on the square of 
the stellar luminosity, SBF magnitudes are especially sen- 
sitive to the brightest stars at a particular wavelength and 
at a given evolutionary phase of a stellar population. As 
shown by various authors SBF magnitudes and, in particu- 
lar, SBF colours, can be used to investigate the properties of 
a specific stellar component in the host st ellar population 
depending on the observing wavelength ('Worthev"l993 
Blakeslcc et al. 2001; Cantiello et al. 2003; Jensen et al_ 
20031: iRaimondo et al.l 120051) . As an example, SBF colours 



involving bluer bands, like B, have been used to study 
the hot stellar comp onent in unresolved systems (e.g., 
ICantiello et al.l l2007b[ ). Likewise, specific phenomena like 
the mass-loss rates in the AGB phase ha ve been analyzed 
by taking advantage of ne ar- IR SBF data (|R aimon dol l200i 
iGonzalez-Lopezlira et al.l l2010t) . As a consequence of the 
quoted relation between the SBF signal and stellar pop- 
ulation properties, we take advantage of the broad pass- 
band coverage to characterize the properties of field stars 
in the galaxy, in order also to provide new constraints on 
the formation history and evolution of the peculiar galaxy 
NGC 1316. 

The organization of the paper is as follows. In iJ2]we present 
a description of the imaging data used, the data reduc- 
tion and calibration procedures. Sections Sj3]and SjDdescribe 
the procedures for SBF measurements and the calibrations 
used to determine the distances. In section [JS]we compare 
the SBF distance to previous results from other indicators. 
We analyze the properties of the unresolved field star pop- 
ulations in ijni and summarize our conclusions in SJT] In 
Appendix El a detailed comparison of SBF and PNLF dis- 
tances is presented. Finally, Appendix IB] presents some de- 
tails on the SBF versus SNe la comparison. 



2 



Cantiello et al.: The distance of NGC 1316 
Table 1. Median ffi and (fh—M) errors from the literature 



Sample 


S{m) 


5(rn-M) Filter 




Number of sources 




(mag) 


(mag) 




(space/ground obs.) 






Optical bands 






ionrv et al. (zUUi) 


U.lo 


A OA T 




2o0 (groundj 


Cantiello et al. f2007a~) 


0.02 


0.09 F814W(- 


'I) 


13 (space) 


Blakeslee et al. (2009) 


0.04 


0.08 F850LP{r 


^ gsDss) 


134 (space) 


Blakeslee et al. f2010) 


0.02 


0.07 F8UW(r^ 


'I) 


9 (space) 






Near-IR bands 






Jensen et al. (1998) 


0.14 


0.19 K' 




16 (ground) 


Liu et al. ('2002") 


0.07 


0.20 Ks 




19 (ground) 


Jensen ct al. (2003) 


0.08 


0.17 F160W (- 


- H) 


79 (space) 



Table 2. Main properties of NGC 1316. 



Alternative names Fornax A, FCC 21, Arp 154, ESQ 357-G 022 

RA(J2000)^ 031i22m41.7s 

Dec(J2000)^ -37dl2m30s 

Galaxy Type^ SO 

Morphological Type^ -1.8±0.7 

Absolute _B-band magnitude^ —22.5 mag 

Recorded SNe la events SN 1980N,SN 1981D, SN 2006dd, SN 2006mr 

cz^ (km/s. Heliocentric) 1788±10 

E{B~Vf 0.021 mag 

iV-I)t 1.132 ± 0.016 mag 

Data retrieved from NED ( |http: / /nedwww.ipac.caltech.edu[ ) ; ^ Hyperleda ( |http://leda.univ-lyonl.fr[ ); ^ ISchlegel et al.l l|1998D : 
* iTonrv et aD (|2001h 



2. Observations and reductions 

This work is based on data of NGC 1316 from the 
Very Large Telescope (VLT) and HST archives. We used 
i) B, V and /-band observations obtained with the 
FORSl Imager at ESO's VLT in Paranal (Program 64.H- 
0624(A), P.L M. Delia Valle), and from the HST ii) ACS 
FA75 W and F85(3LP-b and data from the ACSFCS sur- 
vey (| Jordan et al.l l2007t) . in) WFC3/IR observations in 
the FllOW and F160W filters, plus WFC3/UVIS F336W 
(HST Program ID 11691, P.I. P. Goudfrooij). 

It is useful to note, for the forthcoming discussion, that 
the VLT observations were part of a project aimed at dis- 
covering, monitoring and characterizing the properties of 
the nova population in NGC 1316, with the specific pur- 
pose of deriving the distance to the galaxy using novae 
(jDella V alle & Gilmozzi 2002) . Some relevant properties of 
the target are listed in Table [5J 

In the remainder of this section we describe the re- 
duction and calibration procedures adopted for the data 
from both t elescopes. In all cases we used SExtractor 
(jBertin fc Arn outs 1 996) for the source photometry, and 
the IRAF/S TSDAS task ELLI PSE (based on the method 
described bv lJedrzeiewskillT987[ ) to fit the galaxy isophotes. 

2.1. VLT data 
2.1.1. Data reduction 

The data reducti on was carried out with the VST- Tube 
imaging pipeline (|Grado et al.ll2004 l20li ). specifically de- 
veloped by one of the co-authors of this work (A.G.) for 
data from the VLT Survey Telescope (jCapaccioli et al.l 
120051: ICaoaccioh fc Schipanil [201ll VST). VST-Tube is a 



very versatile software for astronomical data analysis, 
tested against imaging data taken with different tele- 
scopes/detectors, adaptable to existing or future multi- 
CCD cameras (more details will be given in a dedicated 
forthcoming publication, A. Grado et al., in prep.). Further, 
VST-Tube offers the great advantage of fully controlling 
each step of the data processing. 

In Table [3] we report the total exposure times available 
for each VLT filter. Unfortunately, technical problems with 
the camera made a fraction of the total observing time un- 
usable (55%, 60% and 44% of the total exposure time in B, 
V, and I band respectively). The VLT images downloaded 
from the archive and used for this work showed a strong 
degradation during the thirteen observing runs spanning 
nearly two months, from December 1999 to February 2000 
(Figure[T]) . This is partly due to moon illumination, and also 
to the problem of instrument contamination. A decontami- 
nation process of FORSl was p erformed in November 1999, 
and one more in January 2000 (jCavadore et al.lll999D . The 
effect of the last intervention is clearly visible in a ^1.0 
mag jump in the background level between observations be- 
fore and after the decontamination (see Table [7|) . For this 
reason, we decided to select only images taken before 2000 
January 15*'' and, in case of the / filter, we also included 
the last two exposures of 2000 February 4*^* (column "Data 
Used" in Table [7]). The total exposure times used after the 
selection for good frames are also reported in Table [31 

The images were reduced as usual removing the instru- 
mental signatures (overscan correction, bias subtraction, 
fiat field correction and, in the case of the / band, fringe 
pattern removal). The resulting co-added mosaic are ap- 
proximately 6.9 X 7.0 arcmin^, a colour combined image of 
the three mosaics is shown in Figure [2] (upper panel). 
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Table 3. Total and used exposure times for each band. 





B (s) 


y (s) 


/(s) 


Total available 


13200 


9000 


10800 


Used 


6000 


3600 


6000 




Fig. 2. Left panel: false colour combined VLT BVI frames. Middle panel: as upper panel, but for residual frames. The 
location of the ACS frames is outlined in black. Right panel: ACS and WFC3/UVIS combined residual images. Colour 
coding is chosen to enhance the presence of dust. The sites of three over four of the SNe la host by the galaxy are shown. 



2.1.2. Calibration 



The VLT da ta wer e calibrated adopting the 
iGoudfrooii et al.l (|2001bL GOl hereafter) photometry 
as reference, then solving the photometric equations ac- 
cordingly. GOl obtained optical photometry of the sources 
in the field of NGC 1316 using archival NTT/EMMI data 
in 5, V, and / filters. Furthermore, the authors also 
derived JHKs photometry for eight candidate globular 



clusters (GCs) in the galaxy, using the IRAC2 camera 
mounted on the ESO/MPI 2.2m telescope. 

The panels in Figure 12] show the B, V and / VLT 
photometry versus the GOl calibrating data. Note that 
GOl adopted Gal a ctic e xtinction E{B — V) = 0.0 from 
iBurstein fc HeilesI (I1984D while we adopted E{B - V) = 
0.021 from lSchleeel et all (flOOSlfl . 



^ The new meas u remen ts of dust reddening from 
ISchlaflv fc FinkbeineJ l|201ll ) provide ~ 0.002 mag smaller 
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Fig. 3. Comparison between the VLT (this work, t.w.) BVI photometry and NTT data from GOl used for calibrating 
FORSl photometry. 



The comparisons in Figure [3] are obtained in the same 
observational conditions, prior to corrections for Galactic 
extinction. 



2.2. HST data 

HST instruments, thanks to the high resolution and 
the sharp PSF, pr ovide ideal observational datasets for 
SBF analysis (e.g., 'Aih ar et al.l 119971: iJensen et all [200ll 
iBiscardi et all [2008: Blak eslee et al.ll2010[ ). We analvze the 
optical ACS data and provide the first near-IR WFC3/IR 
SBF analysis, with the specific purpose of deriving a con- 
sistent set of ground and space-based SBF measurements 
to secure a reliable distance to NGC 1316. 



2.2.1. NGC 1316 as seen by the ACS 

We analyzed the ACS F475W (- SDSS g band) and 
F850LP (- SDSS z) obser vations of NGC 1316 obtained 
for th e ACSFCS survey fsee lCote et all [200l [Jordan et al.l 
l2007l and references therein for details on the observations 
and the data reductionprocedures for the ACSFCS and its 
twin ACSVCS survey0 

Other ACS observations in the FA35W, F555W, and 
F81AW BV and /, respectively) bands are also available 
from the HST archive. The B-band data are not suitable 
for SBF analysis, while the V and / band have already been 
analyzed by us (Cantiello et al. 2007a), providing results in 
good agreement with the present study (see below). Such 
observations will be used to improve the mapping of the 
dust in NGC 1316. The image processing, including cosmic- 



Galactic extinction, being E{B ~V) — 0.019. The effect of this 
change is at most ~ 0.02, for B magnitudes, and ^ 0.002 mag 
on the colour indices used throughout the present work. 

^ All information on both surveys are also available via web 
at the URL: www.astrosci.ca/users/VCSFCS 



m 0.4 

1 0.2 

f 

i -0.2 
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Fig. 4. Comparison between our photometry and the 
ACSFCS results. Black dots and labels refer to the whole 
sample of ~620 matched sources, gray circles and labels to 
the '~40 selected sources. The median difference and the 
associated rms are also reported. 



ray rejection, alignment, and final image combination, is 
performed with the APSIS ACS data reduction software 
(Bl akeslee et a l. 2003), also used for ACSFCS image align- 
ment. 

For sake of homogeneity with the other datasets used in 
this work, our photometric analysis of the ACS frames was 
performed independently from the ACSFCS one. Hence, it 
is a noteworthy result that our gpirzw and z^ssolp surface 
brightness profiles, derived as described in the following 
section, match within < 0.05 mag the ACSFCS ones from 
Cote et al. (2007). 

The situation is a bit more complex for the photom- 
etry of point-like sources and/or "slightly- resolved" ones. 
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For this comparison we used the preliminary catalog of GC 
ca ndidates from the ACSFCS team, derived as described 
in [Jordan et alj |2004, 2009) , and our photometric catalog 
obtained as described in Section |21 The results of the com- 
parison are shown in Figure SI If the full sample of '^620 
matching sources is used, the median difference between 
our and ACSFCS photometry overlaps with zero within 
the rms scatter. The large scatter between the catalogs is 
due to the independent analysis procedures, especially in 
the way slightly-resolved sources are treated. At the dis- 
tance of the Fornax cluster with the resolution of ACS, 
the GCs hosted by NGC 1316 appear slightly resolved. The 
ACSFCS analysis is optimized to generate accurate pho- 
tometry of GCs with different radii. The aperture correction 
for such sources n eeds to be ev aluated using more refined 
analysis methods (jJordan et a l. 2004, 2009) than the ones 
adopted here (see next section). The gray circles in Figure 
13] show a selection of GC candidates with i) Am < 0.05 
mag and brighter than m = 24 mag in both ACS bands, 
a) galactocentric radius > 60", to reduce the number of 
objects highly contaminated by dust, and Hi) ACSFCS 
estimated radius < 0.03", to select only the most com- 
pact sources, for which the issue of a different treatment 
of the aperture correction should not complicate the com- 
parison. For the selected sample of sources the matching 
is significantly improved (gray symbols in the figure). The 
9F47bW data are statistically consistent in the two cata- 
logs, while a small 0.01 mag offset is seen in the zpsbOLP 
band. However, for the purposes of the present work, such 
an offset only affects the estimate of the contri bution to the 
fluctu ation amplitude due to extern al sources (iTonry et al.l 
[1995 ISode mann fc Thomsen|[T995l : iBlakeslee et al.lll999D . 
Given the level of completeness of the GC catalog, and the 
amplitude of the offset, the impact on zfssolp is negligible. 
As will be shown in section [31 in fact, our SBF measure- 
ments are in very goo d agreement with the ACSFCS ones 
(|Blakeslee et al.ll2009rR . 

2.2.2. NGC 1316 as seen by the WFC3 

The FllOW (- J) and FimW (- H) band images of 
NGC 1316 taken with the WFC3/IR were downloaded from 
the Hubble Legacy Archive together with the WFC3/UVIS 
F336W (~ ?7-band) image. The near-IR images were down- 
loaded for the specific purpose of deriving SBF magnitudes, 
while the Ufssqw has been used to improve the detection 
and masking of dust over the entire set of available images. 

We calibrated the WF C3/IR photometry using the 
VEGA zero points given bv lKalirai et al.l ()2009[ ). To check 
the WFC3/IR photometric data, we made two different 
comparisons, both shown in Figure [5] Figure [5] (left panel) 
plots the WFC3/IR colour Jfiww — Hpwow versus the 
ACS gPAlbW — ZF850LP for the combined ACS-WFC3 data 
set. The crosses give the full sample of sources. In order 
to select the best GC candidates, the black filled circles 
show the sources with zfs5olp < 23.5 mag, photomet- 
ric error Azfs50LP ^ 0.1 mag, SExtractor class-star pa- 
rameter > 0.7 , and galactocentric distance > 45". The 



^ We adopt as reference the photometric VEGA zero points, 
while ACSFCS results use the AB ones. AH ACSFCS data, in- 
cluding the calibration of SBF magn itudes, are transform ed to 
the VEGA magnitude system using ISirianni "eFaLl (|2005l ) zero 
point transformations. 



simple stellar population (SSP) models from the Teramo 
Stellar Population Tools (SPoT, see below) grou iS for the 
age range from 3 to 14 Gyr, and [Fe/H] from —2.3 dex to 
-1-0.4 dex, are also shown in the figure. Even though a non- 
negligible scatter exist between the data of selected GCs 
candidates and SSP models ( ^ 0.2 mag), the match with 
the locus of SSP models is satisfactory for the purposes of 
the present work. 

The right panel of Figure [5] shows the comparison be- 
tween the combined WFC3/IR and GOl near-IR samples. 
The SPoT models are also shown in the figure, and sources 
enshrouded in dust are rejected from the comparison. 
Notwithstanding the large error bars of ground-based obser- 
vations, we find A[{J—H)oqi — {Jfiiqw~Hfi6qw) ~ 0.08] 
mag, in agreement with model predictions. 

3. SBF measurements 

To derive the photometry of sources in each of the selected 
VLT and HST frames, and measure the fluctuation ampli- 
tudes, we used the procedures d escribed in our previous 
works (see ICantiello et al.l l2011aL and references therein) . 
The procedure is basically the same for VLT/FORSl, 
HST/ ACS, and HST/WFC3/IR data with minor differ- 
ences outlined below. 

The main steps of SBF measurement involve: sky back- 
ground determination and subtraction; galaxy model and 
large scale residual subtraction; photometry and masking of 
point-like and extended sources, including dust; power spec- 
trum analysis of the residual frame. We determined the sky 
background by fitting th e surface br ightness profile of the 
galaxy with a Sersic law (jSersidll968h plus a constant term. 
After sky determination, a first model of the galaxy was 
obtained and subtracted from the sky-subtracted frame, 
and a mask of the bright sources was obtained. The large 
scale residuals, still present in the frame after subtracting 
the galaxy model, were removed using the background map 
obtained wit h SExtractor adopt i ng a mesh size ^ 10 t imes 
the FWHM (|Tonrv et al.l IT990t ICantiello et alll200l . In 
the following we refer to the sky, galaxy-model and large 
scale residuals subtracted image as the residual frame. 

The procedure of i) surface brightness analysis and 
sky determination, ii) model fitting and subtraction. Hi) 
sources/dust masking, and iv) large scale residual subtrac- 
tion was iteratively repeated until the residual frame ap- 
peared "flat" in the regions of interest for SBF measure- 
ments, i.e., until the residual did not show any (local) arti- 
fact due to the subtracted galaxy model. The middle panel 
of Figure [5] show the false colour combined BVI image of 
the residual frames. The right panel of the figure, instead, 
shows a combination of the ACS and WFC3/UVIS frames, 
used to map the dust around the center of the galaxy. The 
positions of three of the four SNe la host in the galaxy is 
also indicated in the figure (the region of SN 1980D is not 
covered by either the HST or V LT frames; we adopted the 
revised SNe la coordinates from IStritzinger et al.ll2010l ). 

The surface brightness profiles for all bands analyzed, 
as well as the difference between modeled and observed 
profiles, are shown in Figure [6] 

The photometry of fore/background sources and of GCs 
was derived running SExtractor on the residual frames. 



* See www.oa-teramo.inaf.it/spot 
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Fig. 5. Left panel: colour-colour diagram with the ACS and WFC3/IR magnitudes. Gray crosses mark the full sample 
of common sources, black circles sho w selected GC candidates (see text). Right panel: Comparison between WFC3/IR 
photometry and the IRAC2 data from lGoudfrooii et al.l ()2001a|) . In both panels the predictions from SPoT simple stellar 
population models with different metallicity are shown with different colour/line styles ([Fe/H]=— 2.3, —1.4, —0.7, —0.4, 
0.0 and +0.4 dex shown with dark- yellow, cyan, magenta, green, blue and brown, respectively). The age range is 3 to 14 
Gyr. [See electronic version of the Journal for a colour version of the figure.] 



As described in our previous work, we modified the in- 
put weighting image of SExtractor by adding the galaxy 
model (times a factor between 0.5 and 10, depending 
on the expected ampli t ude of the SBF signa l; for details 
see iJordan et all 120041: iGantiello et~aI1 12005D so that the 
SBFs were not detected as real objects. The aperture cor- 
rection was obtained from a number of isolated point- 
source candidates in the frames an d by making a curve 
of grow th analysis out to large radii (|Cantiello et al.|[2009l 
l2011bD . The outer radius used for FORSl data was 6'.'0. 
For ACS(WFC3) we adopted 0'.'8(1'.'6) and then added an 
extra aperture correction term to infinite radius by us- 
ing t he instrument enci rcled energy tables (for ACS we 
used ISirianni et al1l2005L while for WFC3 we adopted the 
Instrument Handbook, version 4) 

Once the catalog of sources was derived, the next step 
was to fit the luminosity function of the sources, to be 
used to estimate the already mentioned background fluctu- 
ation term due to unmasked faint sources. We obtained the 
fit to the GC and background galaxy luminosity functions 
from the photometric catalog of sources, after removing the 
brightest / saturated point-like sources and the brightest and 
most extended objects. The best fit to the sum of the two 
luminosity functions, and the back ground fluctuation cor - 
rection term, Pr , were derived as in IGantiello et all (l2005l) . 

To measure SBF magnitudes we estimated the az- 
imuthal average of the residual frame power spectrum, 
P{k), then matched it to the power spectrum of a template 
PSF convolved with the mask image, E{k). The total fluctu- 
ation am plitude Pp was obt ained via a robust minimization 
method (jPress et al.lll993) as the multiplicative factor in 



the power spectrum representation P{k) = Pq x E{k) -I- Pi, 
where Pi is the white noise constant term. We used one to 
five different isolated bright point sources in each residual 
frame for the template PSFs. Each PSF, after normaliza- 
tion, was singularly adopted to estimate the SBF signal of 
the galaxy. Finally, the SBF amplitude, P/ = Pq — P^, was 
estimated within circular annuli. The results of the power 
spectrum analysis are summarized in Figure [7l with one 
example for each band. 

The results of the SBF and colour measurements for all 
bands considered are reported in Table 01 For each filter, 
in addition to the statistical error, we report the system- 
atic uncertainty due to PSF fitting. For -ffFieoW: since we 
could only find one good candidate PSF in the frame, we 
assumed a conservative PSF scatter of 0.2 ma ^. In the 
case of Jfiww no good PSF was found over the WFC3/IR 
frame, thus we used a PSF star taken from different ob- 
servations associated with the same HST proposal. As for 
PI^Fi&ow we assumed 0.2 mag PSF uncertainty. 

4. The SBF distance to NGC 1316 

The estimate of distances in the SBF method relies on 
knowledge of the absolute SBF magnitudes. Using the mea- 
surements reported in Table SI together with empirical and 
theoretical calibrations given in Table [5l we obtained the 
distance moduli reported in the (m — M) column of the 
table. In particular we obtained a mean distance modulus 

^ Previo us studie s indica te SBF variations < 0.1 mas with 
PSF (e.g..lLiu et al . 2002; Cantiello et al.ll2005l : iBlakeslee et all 
120101 : IGantiello et al.li2011bl ). 
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Fig. 6. Upper left panel: [/i.336to/WFC3/UVIS, BVI/FORSl, gF475WZF850Lp/ ACS, and JFiiow-ffj^ieow/WFCS/IR 
surface brightness profiles for NGC 1316. Each band is shown with a different line style/colour. A vertical shift is applied 
for sake of clarity as labeled. Lower left panel: As upper left panel, but in i?^/^ scale. The Sersic profi le fitting procedur e 
used leaves the n-index free to vary. Here we use the medi an n = 6 value, which is also consistent with lCote et al.l (|2007[ ). 
The full dots are extrapolated from ICarter et al.l (|l983l ) , accordingly shifted for each band (see colour image) . Right 
panels: the difference between observed and fitted surface brightness profiles. [See electronic version of the Journal for a 
colour version of the figure.] 



of 31.59±0.05 mag and 31.60±0.11 mag with the empirical 
and theoretical equations, respectively. In the following we 
provide some details on the calibration equations adopted. 

4.1. Absolute SBF magnitudes from empirical calibrations 

The empirical calibrations of SBF magnitudes in optical 
bands, in particular in the / band, are the most thor- 



oughly analyze c l (Ffonrv et al.l ll990l 120011: iMei et al.l [20071: 
Blakcsle e et al.l |2010[) . The two aforementioned HST sur- 
veys of the Virgo and Fornax clusters provided an extremely 
accurate calibration of zfsbolp, including some degree of 
non-linearity in the calibration. S ome debate still exists on 
near-IR bands calibrations (see iGonzalez-Lopezlir a et al.l 
2010; and references therein), although, as already men- 
tioned, relevant progres s will be done thanks to the insta l- 
lation of the WFC3/IR (|Blakesleel20'TllFrench et al.ll2012D . 
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Fig. 7. Power spectrum analysis of the BVIzps^qlpJphqwHpiqqiy frames. Each panel shows a different band, as labeled. 
For all bands, the upper panels show the logarithm of the power spectrum of the residual frame (gray dots) and the best 
fit P{k) curve. In the middle panels the difference between observed and fitted power spectra is shown. The flat region 
of log Po{kfit > k) between vertical dashed lines (lower panels) is used to evaluate the best fit parameters Pq, and Pi. 

Table 4. Surface Brightness Fluctuation and colour measurements corrected for galactic extinction. 



VLT data 


(r) 


B-V 


V- 


-I B 


V 


I 


(arcsec) 


(mag) 


(mag) (mag) 


(mag) 


(mag) 


170 


0.84 ±0.03 1.116 ±0.015 33.02 ±0.02 


32.19 ±0.09 


29.81 ±0.07 




PSF uncertainty 




0.09 


0.04 


0.03 


ACS data 


(r) 










z 


(arcsec) 


(mag) 








(mag) 


80 


1.954 ±0.011 








29.11± 0.06 




PSF uncertainty 








0.014 


WFC3/IR data 


(r) 




J- 


-H 


J 


H 


(arcsec) 


(mag) 






(mag) 


(mag) 


80 


1.102 ±0.002 0.72 ±0.01 


27.32± 0.07 


26.28 ±0.04 




PSF uncertainty 






0.2 


0.2 



^ ] Colour data obtained from VLT data using the same masks adopted for WFC3/IR measurements. 



In the upper part of Table [5] we report the distance 
moduli obtained using the empirical calibrations taken from 
literature, together with the adopted calibrations. 

As a first general comment on the empirical equations, 
we must emphasize that the numbers reported in the ta- 
ble are all tied to the same common zero point, i.e., to 
the Cepheid d istances with me t allicity correction to the PL 
relation from \Freedman et al\ \200A) . The resulting zero 
points of the V, I and Hpieow versus V—I calibrations 
are shifted of +0.06, +0.06 and —0.10 mag with respect to 
the calibrations in th e original papers (see appendix A in 
iBlakeslee et"alll2010D . 

For the I vers us B~I calibration, taken from 
ICantiello etaD (|2005h . we do not make any revision since 



the zero point is already based on the chosen set of 
Cepheids. Similar ly, the ^ps.'snLP-band calibration does not 
need any change (jBlakeslee et al.ir2009[ ). 

For the HpiQQw-hand distance estimate, we used two 
i ndepe ndent calibrations. The first derived bv I Jensen et al.l 
(2003) from HST/NICMOS data. We adopted the SPoT 
SSP models to evaluate the changes in the calibration 
due to the difference between the NICM0S2 and the 
WFC3/IR Hpieow passbands. The result is that the 
WFC3/IR Hpiaow zero point is 0.2 mag fainter than the 
NICM0S2 HpiQow one. This is partly expected because 
of the cut at larger/redder wavelengths of the WFC3/IR 
fiher (1400-1700 nm passband, versus 1400-1800 nm for 
NICM0S2). Hence, we add a 0.2 mag to the zero point of 
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Table 5. Distances from SBF measurements. 



Passband 



V 
I 
I 

z 

Jfiiow 



M calibration eq. (m — M) Reference 



Empirical calibrations 
(0.83 ±0.12) + (5.3±0.8)[(1/ 



(-1.68 ± 0.08) + (4.5 ± 0.25)[(1/-/) - 1.15] 
(-1.6 ± 0.1) + (3.0 ± 0.3)[(B-/) - 2.0] 



-2.04+1.41a;+2.60x^+3.722;^ x 



-z) - 1.94 



(-4.8 ± 0.1) + (5.1 ± 0.5)[(y-/) - 1.16] 
-5.17 + 0.70a; + 2.90a;^ x = (g-z) - 1.94 
Weighted Mean 



31.53±0.17 
31.65±0.12 
31.55±0.13 
31.66±0.07 

31.3±0.2 
31.4±0.2 
31.59±0.05 



Theoretical calibrations 



"TTT 

[2] 
[3] 
[4] 

[5] 
[6] 



V 


0.89 + 4.01[(y-/) - 1.15] (0.3) 


31.4 ±0.3 


I 


-1.63 + 5.11[(V-/) - 1.15] (0.3) 


31.6 ±0.3 


I 


-1.63 + 2.65[(B-J) - 2.0] (0.3) 


31.6 ±0.3 


z 


-2.77 + 2.06[(g-2) - 1.94] (0.2) 


31.8 ±0.2 


Jfiiow 


-3.75 + 3.07[(V^-/) - 1.16] (0.3) 


31.3 ±0.3 




-4.86 + 3.59[(V^-/) - 1.16] (0.3) 


31.4 ±0.4 




Weighted Mean 


31.60±0.11 



[Blakcslee et al.' ('2001') 

pTonry et al. (2001) with revised Cepeheid distances 
ICantiello etaP (l2005l) 

Calibration using 'Blakeslee et al.' ('2009'), uncertainties evaluated from eq. (1) in lMei et al.l (|2007f ) 

Calibration from Jensen ct al. (2003) with metallicity correction on Cepheids and NICMOS to WFC3 Hpieow zeropoint 
correction (see text). 

Calibration from lCho. H. et all (|20l3l ). 



the iJensen et al.l (|2003[) empirical relation, assuming a de- 
fault 0.1 mag uncertainty because of the mo del- dependent 
correction term. The second cahbratio n is a preliminary 
result obtained bv lCho. H. et al.l ((2011), based on the ob- 
servations of 16 early- type galaxies in Virgo and Fornax 
specifically obtained to empirically calibrate the SBF for 
the WFC3/IR passband. For the near-IR data, in contrast 
with the optical measurements, the distance modulus in- 
cludes the systematic PSF uncertainty since it is dominant 
with respect to the statistical errors of the SBF measure- 
ment and of the calibration zero point. 

All distances based on the empirical calibrations re- 
ported in Table \b\ agree with each other within the quoted 
uncertainties. The weighted mean of distance moduli is also 
given in the table. 

4.2. Absolute SBF magnitudes from theoretical calibrations 

Various authors have analyzed the possibility to calibrate 
absolute SBF magnitudes using stellar populations synthe- 
sis models, thereby making it a primary d istance indica- 
tor, not l inked to the Cepheids zero point (Buzzoni 1993; 
IWorthevl[T993bl: iBlakeslee et"aLll2001l: [Biscardi et al...200&) . 
In this work we have taken as reference the SBF versus 
colour equations derived using the simple stellar popula- 
tion models from the Teramo SPoT grou p. For a detailed 
review of the SPoT model s we refer to iRaimondo et aD 
(|2005[) an dH aimondol (j2009t ). and references therein. These 
models have been shown to be very effective in match- 
ing the empirical SBF calibration in different bands, as 
well as in reproducing the resolved (colour magnitude di- 
agrams) and unresolv ed (colours, rnagnitudes) properties 
of stellar populations (jBrocato et al.l[2000t ICantiello et al.l 



l2007al:ICantieIioll2012l) . We used the updated version of the 
SPoT models, which for the photometric bands and chem- 
ical composition used in t his section confirms the results 
obtained from the previous IRaimondo et al.l (j2005l) models 
(G. Raimondo, private communication). The grid of mod- 
els used has [Fe/H]= —0.4, 0.0, +0.4 dex and ages from 3 
to 14 Gyr. The choice was made based on the age and 
chemical composition properties of fields stars in NGC 1316 
derived from various independent spect r o-photometric in- 



dicators (e.g. iTerlevich fc ForbesI [200l ISilva et al.l [20081: 

Konami et al. 20101), and also confirmed by GC analysis 
(Goudfrooii et al.l 200 1a.b. sec also Under these as- 
sumptions, and using a bootstrap approach, we obtained 
the calibration equations reported in the lower part of Table 
[5] For each equation, the scatter of the M versus colour re- 
lation is also tabulated. The uncertainties on distance mod- 
uli are derived by summing in quadrature the scatter of the 
theoretical calibration and the uncertainty on ifi. 

The distance moduli derived with the theoretical cali- 
brations from the SPoT models are given in the TablejS] All 
reported distances agree to within the quoted uncertainties. 

4.3. Combining SBF-based distance moduli (plus a 
fundamental note on uncertainties) 

The SBF distance of NGC 1316 presented in gUand 
is based on a self-consistent treatment of empirical and the- 
oretical SBF calibrations. Empirical calibrations have been 
tied to the same common zero point reference. On the other 
hand, the theoretical calibrations are based on a well de- 
fined set of SSP models ranging from V to i?Fi60W, i-e., 
a wavelength interval that, in terms of SBF and colours, 
samples very different stellar population (sub)components. 
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Fig. 8. V, I and i?_Fi6ow-band SBF versus {V~I) pre- 
dictions from five independent stellar population synthe- 
sis codes: solid line (black), dotted line (red), dashed 
line (blue), long-dashed line (gr een), dot-dashed (cyan ) 
line show the u pdated SPoT, [B lakeslee et al." (2001'), 



Liu et al.l (I2002D. iMarin-Franch fc Aparicio (10Q6), and 
Gonzalez-Lopezlira et al. best fit lines, respectively. 



Empty (magenta) circles show the corresponding empirical 
relations, also given in Table [5] [See electronic version of 
the Journal for a colour version of the figure.] 



sphere models, etc.), and then analyze the effects on SBF. 
A rough (indirect) estimation of the uncertainty on the- 
oretical calibrations is to compare the results from inde- 
pendent models, obtained from different SSP codes relying 
on independent physical input and algorit hms. A first at- 
tempt along these lines was carried out by ICantiello et al.l 
(|2003f ). To estimate this uncertainty for the present work, 
we com pare SSP model predictions from Blak eslc e et al.l 
(2001), Liu et al.l (I2002D iMarfn-Franch fc Apariciol (|2006l ). 
'Gonzalc z-Lopezlira et al.n 20101 ) and the SPoT models. For 
all these models, we obtain linear fits to the V, I and 
Hpieow SBF amplitudes versus {V—I) colour (Figure [5]). 
The empirical calibration for each band is also shown in 
the figure. Note the different range of colour used for each 
equation, depending on the range of vali dity of the em- 
pirical equation ( Tonry et al. 2001; Blake slee et al.l 120011 : 
I Jensen et aLll2003t ICantiello et al.ll2007aD . We find that a 
scatter of ~0.2 mag provides a first approximation to the 
systematic uncertainty of SBF calibrations derived from 
SSP models. A discussion on the origin of the scatter be- 
tween models is beyond the scope of this paper. We refer 
the interested reader to the quoted papers and references 
thereafter. 

In conclusion, by combining the weighted mean distance 
moduh, {fh — M) empirical = 31.59 ± 0.05(stat.) ± 0.20(sys.) 
mag and {m-M)theoreticai = 31.60±0.11(stat.)±0.20(sys.) 
mag, we conclude that our best estimate of the distance of 
NGC 1316 is (m - M) ^ 31.59 ± 0.05(stat.) ± 0.14(sys.) 
mag, or d = 20.8 ± 0.5(stat.) ± 1.5(sys.) Mpc. 



It is noteworthy that the empirical and theoretical evalu- 
ations agree very well with each other. Taking into account 
that the two methods are based on independent calibration 
procedures, subject to different types of systematic and sta- 
tistical uncertainties, this result suggests that both types of 
errors are reasonably well constrained. 

For what concerns the systematic errors in our measure- 
ments, summing up all expected sources of uncertainty for 
each one of the three instruments considered - filter zero 
point, data reduction, calibration zero point, PSF normal- 
ization - the expected systematic uncertainty is ^0.1 mag, 
with the exception of near-IR bands, where the contribu- 
tion from the PSF normalization is dominant with respect 
to the others. 

For the systematic errors in the empirical calibrations, 
recall that all relations are li nked to the same Cepheid zero 
point (jFreedman et al.ll2001h . which accounts for a further 
systematic ^ 0.2 mag uncertaintj0- 

Estimating the systematic uncertainty in theoretical 
SBF calibrations is not a simple task. One way to get 
such an estimate would be to change the ingredients in 
the SSP code (stellar tracks, initial mass function, atmo- 



A further component to the systematic er ror comes from 
bandpass mismatch with various telescopes. iBlakeslee et al.l 
l|2001l ) presented a discussion of this issue (see their Sect. 5.5, 
Fig. 15 in particular), showing that at the colour of the SBF, 
(V-I) - 2.4 for NGC 1316, the difference in the standard 
Cousins / and the HST / can vary by ±0.02 mag. 



5. Comparison with distances from the literature 

As already mentioned, the VLT observations used here were 
obtained to detect and study the population of novae in 
NGC 1316, and to use them t o derive the galaxy distance. 
iDella Valle fc Gilmozzil (pOOl . using the data of four de- 
tected novae (the first detected beyond the Virgo cluster 
at the time), and the Buscombe-de Vaucouleurs relation 
(B uscombe fc de Vaucouleur£lll955l : ICapaccioli et al.iri989l 
I1990l) " set an upper limit on the distance of the galaxy equal 
to 22.4 Mpc, and a lower limit of 16 Mpc from the data 
of Nova A in their sample. In spite of the difficulties (ob- 
servational and statistical) in using the novae in Fornax 
A to estimate the gala xy distance, the range obtained by 
IDella Valle fc Gilmozzil is in excellent accord with our esti- 
mates. Moreover, the lower limit distance is also in complete 
accord with the SNe la as well as the PNLF distances (see 
below). 

In the following sections we compare our distance esti- 
mates with others available in the literature. We consider 
the results obtained from three different distance indica- 
tors: Type la supernovae, the planetary nebula luminosity 
function (PNLF), and mean properties of the GC system. 
Other distances indicators, such as the TuUy-Fisher rela- 
tion or the fundamental plane, cannot be used reliably for 
NGC 1316 because of its irregular post-merger morphology; 
moreover, these indicators are normally applied to groups 
or clusters as a whole, rather than giving precise distances 
to individual galaxies. 
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5.1. Comparison with SNe la 



lAihar et alJ (|200lf ). and later ICantiello et alJ (|2011a[ ). have 
presented a comparison of SBF and Type la SNe distances 
for a total of 15 different SNe la in 14 galaxies. Both stud- 
ies found excellent overall agreement between the two dis- 
tance indicators, including the case of NGC 1316, provided 
that a co nsistent set of Cephei d-bas ed distances is used. 
Simila rly, iFreedman e t al. (200l|) and iFreedman &: Madord 
(|201C() derive essentially identical values of Hq from these 
two methods when calibrated consistently via Cepheids. 

NGC 1316 has been a prolific producer of Type la su- 
pernovae, with four recorded events: SN1980N, SN1981D, 
SN2006dd, and the fast declining SN2006mr. Ah SNIa 
light-curves have been used to determine the distance of 
the host galaxy using several methods. Furthermore, being 
one of the nearest bright galaxies with well sampled SNe la 
light-curves, NGC 1316 is frequently use d in calibrating 
samples for high- redshift SNe la (e.g., iJha et all l2007t 



iBurns et all 1201 ll). B ased on the MLCS SNIa distance 
method (jRiess et all I1998D . iGoudfrooii et all (|2001aD re- 
ported (m-M) = 31.80 ±0.05 mag for NGC 1316 and con- 
cluded that it was therefore ~ 0.25 mag m ore distant than 
the re st of the Fornax cluster, for which iFerrarese et al.l 
(|2000tl gave a mean distance modulus of 31.54 mag from 
Cepheids and other indicators. However, no detai l s on th e 
calibration, etc., w ere given by Goudfrooii et al.l (|2001a^ . 
lAihar et all (l200lD re p orted both MLCS and Amis dis- 
tances (iPhiUipsI 119931: iHamuv et all MM) for SN 1980N 
in NGC 1316 under different calibrations. Rescaling their 
results to Hq = 72 km Mpc~^ gives (m—M) — 
31.30 ± 0.15 and (m-M) = 31.51 ± 0.10 from the MLCS 
and Amis methods, respectively. Interestingly, the Amis 
value coincides ex actly with the mean SB F distance for the 
Fornax cluster by'Blakeslee et al.' (20091). 

On the other hand, Stritzinge r et all ()2010l StrlO here- 
after) have recently reanalyzed the SNIa distance to 
NGC 1316, and obtained (m - M) = 31.25 ± 0.03 (stat.) 
±0.04 (sys.) mag, based on the analysis of the three "nor- 
mal" SNe la (adopting Hq = 72 km s"^ Mpc"\ see Table 
[6]). The authors also obtain values of (m — M) ~ 31.7 — 31.8 
mag from the data of the fast declining SN2006mr; how- 
ever, they consider as doubtful the suitability of fast- 
declining SNe la for estimating distances. Thus, the best 
distance modulus to NGC 1316 obtained by StrlO is ^0.34 
mag fainter than ours, equivalent to a 17% smaller distance. 
Note that the new distances derived by StrlO, as pointed 
out by the authors, are based on a thorough analysis and 
discussion of the four SNe la, while little or no detail were 
given in previous studies. Hence, such a difference is very 
interesting, and needs further analysis, especially because 
of the primary role in the cosmological distance scale of 
both these indicators. 

As a first general comment, it is useful to highlight that 
two of the normal SNe (SN 1980N and SN 1981D) were ob- 
served in the pre-CCD era. So, although they afford a nice 
opportunity to test the internal consistency of SNe la dis- 
tances, one should not forget that the quality of the data 
of these objects is lower than what is currently being ob- 
tained. Bearing this in mind, we proceed to our analysis 
taking the SNe la distances, and their associated system- 
atic/statistical errors, from the cited references. 

The SNe la distances to NGC 1316 obtained before the 
study of StrlO, are listed in Table [HI For consistency with 



our and StrlO's measurements, all distances are converted 
to a scale with Hq = 72 km s^^ Mpc~^. In the table we 
do not report the systematic errors, which are not given 
by all authors. Taking the weighted mean of all measure- 
ments in Table m except those of StrlO, adopting 0.4 mag 
default uncertainty where no error is reported, we obtain 
(m — M) = 31.35 ± 0.05 mag. Although such value agrees 
with our estimates better than the StrlO result, we caution 
the reader against these "general" averages, reported here 
only to emphasize the relative distribution of distances with 
respect to a reference point. Some of the values reported 
in Table [51 in fact, are obtained using the same meth- 
ods/objects but under different assumptions (e.g. MLCS 
and MLCS2k2), and/or with different calibrators, so that 
the reported average does not necessarily have a correct 
physical meaning. 

StrlO derives the distance to NGC 1316 using three dif- 
ferent methods: the EBV, the Tripp method, and the near- 
IR light-curves of Type la SNe. Each method has its own 
calibration, and is quantitatively independent from the oth- 
ers. Hence, there can be various possible causes for the dif- 
ference between our and StrlO distances. First, for the EBV 
method, the authors adopt negligible internal extinction for 
all four SNe la, but also comment that the spectroscopic 
analysis provides results "totally inconsistent with the low 
host galaxy reddening" (StrlO, Section 4). However, the op- 
tical and optical/near-IR colours of these objects are con- 
sistent with minimal to unreddened supernova (Figure 13 
of StrlO). Second, for the near-IR method, which is intrinsi- 
cally much less affected by the ho st internal extinct i on, w e 
find that the calibration used, from lKrisciunas et al.l (|2009[ ). 
is based in part on a compilation of Cepheid and/or SBF 
distances that is not internally consistent (see the discus- 
sion in Appendix |B| . 

Interestingly, StrlO obtain a SN la distance modulus 
for the Fornax cluster member NGC 1380 of (m—M) = 
31.611 ±0.008 mag, which is consistent wit h the SBF result 
for th e same galaxy of 31.632 ± 0.075 mag ()Blakeslee et al.l 
I2009D . and very similar to the SBF distance for NGC 1316. 
In fact, there is little significant variation in the SBF dis- 
tances among the magnitude-limited sample of 43 early- 
type Fornax galaxies studied bv Blakcslee et al.l (j2009D . and 
most are consistent with the mean Fornax SBF modulus of 
31.51 ±0.03 ma£|. 

Finally, we must note that the iTripd . the maximum 
near-IR magnitudes, and the MLCS2k2 methods applied 
to the fast declining SN2006mr give distance estimates in 
agreement with the ones presented here. In spite of this, 
it should also be emphasized that the debate on whether 
fast declining SNe la can be used to determine precise dis- 
tances is still open (so that different autho rs include or not 
these objects in thei r final saniples; e.g. iJha et all 120071 : 
.FolateUi et al.,,2010; .Bmns et aLlboiU StrlO). 



T he recent Type la supernova SN2012fr (|Childress et al.l 
120121 ) occurred in NGC 1365, a giant barred spiral galaxy 
in the direction o f the Fornax cluster w i th a measured 
Cepheid distance (|Silbermann et al] 119991 : IFreed man ct al.l 
200 li ). However, from a comparison of Cepheid and SNc la dis- 
tances, iSuntzeff et all fl999) suggested that NGC 1365 IS ac- 
tually ~ 0.3 mag in front of the Fornax cluster eUipticals; 
iKelson et al"] (|2000l ) came to a similar conclusion based on the 
fundamental plane. 
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Table 6. Type la Supernova distances to NGC 1316 available in the literature. 



SN 


(m-M) 




Method 


References/Notes 




Measurements from StrlO 




oUlN+oiJJH-UDaa 


31.180 ± 0.013 




SNooPy EBV 


Author's best estimate 


80N+81D+0Ddd 


31.25 ± 0.03 




Tripp 


Author's best estimate 


oUiN -l~olU-rUDQa 


31.203 ± 0.012 


(m 


- M)^a:, NearTR 


Author's best estimate 


bIN ^UUbmr 


31.83 ± 0.07 




Tripp 


Flagged as doubtful 


oiM ZUUDllil 


31.739 ± 0.005 


(m 


- M)^a^ NearTR 


Flagged as doubtful 




Other Measurements^ 




SN 1980N 
SN1981D 


31.45 
31.35 






Hamuv et al. (1991) 
Hamuv et al. (1991) 


SN 1980N 


31.44 ±0.32 




MLCS 


Kriscmnas et al. (2000) 


SN 1 981 n 


31.03 ±0.3 




MLCS 


Krisriiinat; et al r^OOO"! 


SN 1980N 


31.30 ±0.15 




MLCS 


Aihar et al. (2001) 


SN 1980N 


31.51 ±0.10 




Amis 


Aihar et al. (2001) 


SN 1980N 


31.35 




■^max 


Roindl et al. (2005) 


SN1981D 


30.98 




^^max 


Roindl et al. (2005) 


SN 1980N 


31.38 ±0.13 




MLCS2k2 


Jha ct al. (2007) 


SN1981D 


30.78 ±0.24 




MLCS2k2 


J[ha ct al. (2007) 


SN 2006dd 


31.34 ±0.11 




MLCS 


Maoz & Mannucci (2008) 



All distances are converted to a scale with Hq = 72 km s ^ Mpc 



5.2. Distances from Globular Cluster System properties 

Distances to NGC 1316 derived from the properties of the 
globular cluster system have generally been flagged as un- 
reliable by the authors due to the peculiar properties of the 
galaxy and its GC system. 

iGomez et al.l (|200l[ ) made one of the first attempts 
to constrain the distance to NGC 131 6 with the G lobular 
Cluster Luminosity Function (GCLF, lHarrisl[200l[ ). using 
ES0/EF0SC2 data. The BVI weighted average distance 
modulus they provided is (m — M) ~ 31.4 mag, based 
on a calibration consistent with ours. However, the au- 
thors specifically commented on the existence of red and 
blue sub-populations of GCs and concluded that their cat- 
alog was not sufficiently deep to reliably constrain the 
galaxy distance. In an earlier study of five bright Fornax 
cluster galaxies, ? had noted that overall the GCLF of 
NGC 1316 was not well fitted by a Gaussian (reduced 
Xn > 3, as compared to Xn ~ 1 for the four others), 
and thus could not prov ide a rehable GCLF distance for it. 
IGoudfrooii et al.l (|2004f ) presented a deeper analysis of the 
GC system based on ACS data. These authors confirmed 
the previous results of a blue GC sub-population, consistent 
with a Gaussian luminosity function, and a further com- 
ponent of red GCs with a power la w luminosity function. 
More recently, iVillegas et all (l2010l) for the ACSFCS sur- 
vey, write "NGC 1316 is [...] not included in the fits because 
the observed GC system in this galaxy is highly influenced 
by its interaction and proximity with its satellite galaxies, 
and therefore our GCLF fit is not reliable. " 

Masters et "all (|2010[ ) use the GC-radii method 



([Jordan et al. 20051 ) to analyze the distances of Fornax 



Cluster galaxies. However, the authors warn that, while 
the method is effective for typical GC systems, it cannot be 
applied to NGC 1316 due to the large number of extended 
GCs. 

In conclusion, the GC-based distances reported, for ex- 
ample in the NED archive, (m — M)GCLF,g — 33.55 mag, 
(m - M)gclf,z = 33.68 mag, and (m - M)Gc-radius = 
30.59 ± 0.11 mag, are excluded from the comparison. 



5.3. Planetary Nebulae Luminosity Functions 

iFeldmeier et al.l (|2007t ) reported a PNLF distance mod- 



ulus to NGC 1316 of 



M) 



31.26 



-1-0.09 
-n 19 



lag. 



This 



result relies on the PNLF calibration bv ICiarduUo et al 
(2002), whose zero point is tied to the IFreedman et al 



(|2001i ) sample of Cepheid's distance moduli with no de- 
pendence on metallicity of the PL relations. As discussed 
above ( H4.ip . the SBF calibrations presented here are based 
on the metallicity-corrected Cepheid PL relations. Hence, 
to properly compare the PNLF and SBF distance, a c orrec- 
tion is necessary. In a recent review, iCiardullol (|20T1 finds 
that if one adopts the metallicity-dependent Cepeheid PL 
relations, the PNLF zero point brightens by 0.07 mag, going 
from M* = -4.46 ±0.05 to M* = -4.53 ±0 . 04 ma g (exter- 
nal scatter a = 0.16 mag). IFeldmeier et al.l (|2007l) actually 
used M* — -4.47; thus their PNLF distance modulus for 
NGC 1316 becomes (m-M) = 31.321°;'^^ ma^. 

Both PNLF distances given above agree within quoted 
uncertainties with the StrlO distance. However, the SBF 
and the updated PNLF distances are consistent within 
the given statistical and systematic errors, notwithstanding 
that the difference between them remains non-negligible. In 
other words, within the given uncertainties, the SBF dis- 
tance is consistent with the PNLF, and the latter with the 
SNe la, but SBF and SNe la are not consistent with each 
other. Moreover, with the revised zero point that corrects 
for the metallicity dependence of the Cepheids, the other 
PNLF distances published in the past decade for Fornax 
galaxies are (m-M) = 31.10tE! }^ mag for NGC 1380 
(jFeldmeier et al.ll2007l) and (m-M) = 31.46 ± 0.18 mag 
for NGC 1344 (jTeodorescu et al.ll2005D . Thus, although the 
PNLF method finds systematically lower distances to both 
the Virgo and Fornax clusters, it is consistent with SBF 
in placing NGC 1316 well within the distance range of the 
other Fornax cluster galaxies. Further details on the SBF- 
PNLF comparison are given in Appendix |21 



We have not taken into account a further —0.04 mag cor- 
rectio n term to SBF distance moduli cited by ICiarduUo et al.l 
(|2002l . their section 6). Including such an extra correction to 
our distance modulus implies {fh—M) ~ 31.55 mag. 
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6. SBF and integrated colours to constrain stellar 
population properties 

The comparison of model predictions with observed SBF 
magnitudes and SBF colours to understand the prop- 
erties of the host galax y has already been successfully 
used by different authors ("To nrv et alJ 1990:'Buzzon illl993l: 
Jensen et al. 2003; Raimond o et al.l [ 2 005: Cantiell o et alJ 
2007bt IBuzzoni fc Gonzalez-L6Dezlirall200a) . 

The SBF colour versus integrated B-V for NGC 1316, 
plus two more Virg o cluster galaxies, N GC 4374 and 
NGC 4621 (data from ICantiello et al.ll2011aD . are shown in 
Figures [9] and [TOl The data in the figures are compared to 
SPoT SSP models computed with standard assumptions for 
different metallicity and age (panels a — c), while the pan- 
els from d to n show the predictions for solar-metallicity 
standard and non-standard models: specifically, i) in pan- 
els d — f SSPs with an enhanced hot horizontal branch 
component (HHB, having ^ 50% stars in the canonical HB 
and ~ 50% HHB stars) are considered; ii) in panels g — i 
results obtained by adding a fraction of very young stars of 
30 or 100 Myr to the old, solar [Fe/H] component (in pro- 
portion 1:1000) are plotted; and Hi) in panels from I to n 
predictions for models with a mix of an old solar metallic- 
ity component, and a further one with [Fe/H]=— 2.3, —1.3, 
—0.7, —0.4 and -1-0.4 dex (in portions 1:1) are shown. 

Panel c in Figure clearly shows the well-known 
age-metallici ty degeneracy that affects classical integrated 
colours (e.g.. lWorthevll99^ . The position of galaxies in this 
panel overlaps nicely with SSP models, though nothing can 
be said about the stellar content. Because of the overlap 
between models, the field stellar component in NGC 1316 
could be either older and more metal poor or younger and 
more metal rich than the other two galaxies. In contrast, 
the SBF colour versus B—V models shown in panels a — b 
are rnuch less affected by the degeneracy, especially for the 

V — Hpieow colour. The data for the three galaxies lie near 
the region of solar metallicity models. A possible interpre- 
tjition of the relatjve positions of the three galaxies in the 

V — HpiQQw and V — I versus {B—V) colour planes, is that 
NGC 1316 hosts a field component that is as metal rich as in 
NGC 4621, with [Fe/H]~0.0 dex, but significantly younger 
{3 < t (Gyr) < 10). In fact, while both the galaxies in 
Virgo lie near the edge of old SSP models, this is not the 
case for NGC 1316, for which the stellar light is polluted by 
an intermediate-age component. This is seen more clearly 
in some of the other panels of the figure. 

Inspecting panels I — n we find that, while the colour- 
colour panel n shows the expected age-metallicity degener- 
acy, the situation changes for SBF colours (panels I — m). 
The Virgo cluster galaxies, in fact, overlap with the region 
of models obtained with the mixing of old t=14 Gyr SSPs 
with different [Fe/H]s. In both the I and m panels NGC 1316 
lies above the line of mixed old SSPs, suggesting that a 
younger SSP is necessary to obtain a good match with 
the models. Again, this is not surprising since NGC 1316 is 
a known example of an intermediate- age merger remnant, 
and also has Mg2 ^ 0.25 mag, while the other two targets 
have Mg2 > 0.28 mag. 

The data to models comparison is less straightforward 
when i3-band SBF magnitudes are consider ed. As discussed 
by oth er authors (e.g., Worthey 1993^ ICantiello et al.l 
l2007bf) . B-band SBF cannot be used to get reliable galaxy 
distances, both because the amplitude of the signal is 2-3 



orders of magnitude fainter than in optical/near-IR bands, 
and because of the strong sensitivity to stellar population 
properties. This is depicted in Figure [TU] (panels a — c), 
where we plot SBF colours obtained with B versus the 
B—V for the three galaxies, and the standard SPoT mod- 
els. One major difference with the results in Figure [S] is 
the substantial mismatch between data and models seen in 
panel b and, especially, in panel c. The B—V versus B—V 
models are clearly affected by the age-metallicity degener- 
acy, but in this case the data are offset by ^ 0.6 mag with 
respect to the models. The mismatch is less evident for the 
B—HpiQQw colour, due to the much larger baseline of this 
colour. 

If non-standard SPoT SSP models are taken into ac- 
count (panels from c? to n in Figure [TU)) we notice that: 

— models with an enhanced number of hot HB stars pro- 
vide a good match to the data; 

— the mean metallicity of the dominant stellar component 
in NGC 1316 is apparently lower than that of the two 
Virgo galaxies. The models shown in panels d — /, in 
fact, are obtained assuming HHB "enhancement" for 
three different ages and solar metallicity. Hence, taking 
into account SSP models with lower [Fe/H], i.e., bluer 
colours, the matching with the position of NGC 1316 
will improve; 

— the presence of a diffuse very young stellar component, 
with an age of 30 Myr, seems to reduce the mismatch 
with data (panels g — i). However, even in the case of the 
merger remnant NGC 1316, the presence of such a young 
diffuse stellar component is unlikely even if the pres- 
ence of a relatively young stellar population is expected 
in strong radio-emitter galaxies ([Delia Valle fc Panagial 
2003); 

— the improved data and models matching in panel I seems 
to support the possibility of an old [Fe/H] mixed stellar 
component. However, this is ruled out by the compar- 
ison in panel n, showing a poor match to the B—V 
predictions for the same data set. 

At this stage it is useful to recall that the HHB scenario 
is supported by other independent obser vati ons. First, we 
mention the results by Bro wn et al.l ([200 0^) and'B rown et al.l 
(2008) who found a significant fraction of HHB stars in 
M32 using HST UV data. Second, the puzzling pres- 
ence of a strong UV emission in some regular early-type 
galaxies, discovered several decades ago ([Code et al.lll972l : 
Ber tola et allll980|) . is now widely interpreted as the pres- 
ence of an old hot stellar component. Although the mech- 
anisms regulating this component are not well un derstood 
(iPark fc Ledll997l: iKavirai et ai1l2007l: [Halll2008ll ,_some of 
these old hot stellar sources may have effects on B, as in- 
depen d ently predicted by v arious SSP models ([WorthevI 
Il993al : [ Cantiello et al.l 20031 ) , or based on empirical evi- 
dence (IShopbell et al. [ [iggalSodemann fc ThonlsCTi [Tggi 
ICantiello et al.ll2007bD . 

We emphasize that, when computing the SBF calibra^ 
tion equations using the models with HHB, we find that B 
brightens by ^^0.2 mag at fixed colour. The effect is ^ 0.05 
mag in other bands (negligible redwards of z). Such be- 
havior, again, highlights the uselessness of B for distance 
determinations, and confirms the interest on SBF in blue 
bands for analyzing the properties of unresolved blue hot 
stellar components. 
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Fig. 9. Panels a — c: SPoT standard SSP models compared with SBF measurements of NGC 1 316 (full black circle). The 
data of NGC 4374 (filled square), and NGC 4621 (filled triangle) from ICantiello et alj (|2011a|) are also shown. Different 
line styles mark different [Fe/H] contents: dot-dashed (magenta), long-dashed (green), solid (blue) and dashed (drown) 
refer to [Fe/H] =—0.7, — 0.4, 0.0 and 4-0.4 dex, respectively. Models in the age-range 3-14 Gyr are shown. The arrow in 
panel a indicates the direction of increasing age. Panels d ~ n: SPoT models obtained with non-canonical assumptions. 
In all cases the initial (reference) population has solar metallicity, t~ 14 Gyr, and is shown with a blue star, while the 
final composite population is shown with a different colour and connected with a dotted line to the reference model. 
The symbols for observational data are the same as in panels a — c. Panels d — f: A fraction equal to 50% of the total 
HB-stars is simulated being HHB. For these models three ages are considered (t=10, 12 and 14 Gyr, increasing ages are 
marked with larger symbols). Panels g — i: A young population of t=30 Myr (red), or t=100 Myr (magenta) is added to 
the old solar one. The fraction in mass of old to young stars is reported in the lower panel. Panels / — n: old SSPs with 
various [Fe/H] are mixed to the solar one, as labeled. The mass fraction metal-poor to standard SSPs is shown in the 
lower panel. [See electronic version of the Journal for a colour version of the figure.] 



In conclusion, the present analysis of the stellar popula- 
tion properties for NGC 1316 seems to confirm the results 
of NGC 4374 and NGC 4621, i.e., that a diffuse component 
of hot old stars contributes to the SBF signal in the B- 
band. Furthermore, the relative comparison of SBF and 
colour data of NGC 1316 with those of the two galaxies in 
Virgo seems to indicate that the dominant stellar compo- 
nent in NGC 1316 is younger and slightly less metal rich, 
as expected for an intermediate-age merger remnant. 



7. Summary 

We have measured SBF magnitudes in NGC 1316, the 
brightest galaxy in the Fornax cluster, using ground-based 
VLT/FORSl data in BVI-hands, and space-based gF475W 
and 2_F850LP from ACS/WFC, plus Jpiiow and i?Fi60W 
WFC3/IR observations. 

The distance of NGC 1316 is particularly interesting in 
the context of the cosmological distance scale. The Fornax 
cluster, in fact, is the second largest cluster of galaxies 
within i$ 25 Mpc after the Virgo cluster. However, in con- 
trast to Virgo, the line-of-sight depth of Fornax is small, 
enabling accurate calibration of distances without the ad- 
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Fig. 10. Same as in FigurelHl but B is considered instead of V. [See electronic version of the Journal for a colour version 
of the figure.] 



ditional scatter intrinsic to the spatial extent of the clus- 
ter. Furthermore, NGC 1316 is among the galaxies with the 
largest number of detected Type la supernovas (SN 1980N, 
SN1981D, SN2006dd and SN2006mr). For this reason, it 
is a unique place to test the consistency of SNe la distances, 
both internally and against other distance indicators. 

Using our SBF measurements in VIzfs5olpHfi60W 
and available empirical calibration of the absolute M, we 
obtained a weighted mean distance modulus to NGC 1316 
(m - M)empiricai = 31.59 ± 0.05(stat.) ± 0.20(sys.) mag. 

Additionally, we obtained SBF distances from 
VIzFa5QLpJFiiowHFiGQw data based on theoreti- 
cal calibrations derived from the SPoT SSP models. 
The resulting weighted mean distance modulus is 
(m - M)theoreUcai = 31.60 ± 0.11(stat.) ± 0.20(sys.) mag. 

The good agreement between SBF distances obtained 
from empirical and theoretical calibrations is a notable re- 
sult. The two classes of calibration are, in fact, completely 
independent of each other, one relying on the first two rungs 
of the cosmic distance scale, the other on the present knowl- 
edge of the various ingredients that go into stellar popula- 
tion synthesis (stellar evolution theory, stellar atmospheres, 
etc.). Furthermore, since SBF magnitudes over such a wide 



range of wavelengths depend on the properties of stars in 
different evolutionary stages, even the agreement between 
distance moduli obtained in different bands should be re- 
garded as a remarkable result particularly in terms of the 
high degree of reliability reached by SSP models. 

By combining the distance moduli from the two types 
of calibrations, we obtain (m — M) = 31.59 ± 0.05(stat.) ± 
0.14(sys.) mag, or d = 20.8 ± 0.5(stat.) ± 1.5(sys.) Mpc. 

This distance modulus agrees generally well with es- 
timates obtained from other indicators, and with SNe la 
light-curves analysis obtained before 2010. A non-negligible 
difference exists with the m ost recent analysis based on 
SNe la by IStritzinger et al.1 ([2010,1 , who obtained a best 
estimate for the distance ~ 17% smaller than ours. The 
possible sources of the disagreement may be related to the 
complex issue of the internal extinction, and to zero point 
calibration issues of both distance indicators. 

The comparison to PNLF is also subject to the linger- 
ing problem of zero point calibration. When placed in a 
Cepheid distance scale consistent with ours, the PNLF dis- 
tance to NGC 1316 is d = 18.4 ± 1.0(stat.) ± 1.5(sys.) Mpc 
and agrees within the quoted errors with our SBF distance. 
However, the difference between SBF and (the updated) 
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PNLF distance remains non-negligible, and lowering the 
uncertainties (especially systematic; see Appendix |2| 
would be a desirable result for both indicators. 

In order to analyze the properties of the dominant stel- 
lar component in the galaxy, we compared SBF colours and 
integrated colours to SSP model predictions. We found that 
the stellar light of the galaxy seems to be dominated by a 
[Fe/H] ^ 0.0 dex, intermediate age stellar population. The 
comparison with analogous measurements for NGC 4374 
and NGC 4621 also supports a scenario in which field stars 
in NGC 1316 have a younger age and slightly lower metal- 
licity than the two bright Virgo cluster members. 

Moreover, we found that SBF predictions from standard 
SSP models do not match with observations if _B-band 
SBF data are included in the model comparison. As in 
the cases of NGC 4374 and NGC 4612, which showed a 
similar mismatch to models, we used the SPoT stellar 
population synthesis code to generate SSP models with 
non-canonical properties. In particular, we considered the 
following three cases: starting from an old t ~ 14 Gyr 
population with solar metallicity we have 1) enhanced the 
content of hot HB stars, 2) added a very young diffuse 
secondary component, and 3) added a more metal poor 
SSP. As in the previous case (jCantiello et al.ll2011a[l . the 
simulations seem to favor the HHB component scenario. 
Assuming a contribution to B from such hot HB stars 
removes the discrepancy between the data and models in 
this band, yet has negligible effect on SBF in other bands, 
i.e., it does not affect the theoretical calibrations used to 
obtain distances. 

Our results on the distance and stellar population prop- 
erties of NGC 1316 based on SBF analysis have shown that, 
despite the great progress in recent years, many issues re- 
main open on both topics. Concerning distances, the cali- 
bration of distance indicators, and the treatment of error 
propagation in the distance scale, still need to be accurately 
and consistently analyzed. Concerning stellar populations, 
SBF colours, as a new and independent stellar population 
analysis technique, seem to provide useful constraints to 
the properties of field stars, hidden to many classical pho- 
tometric indicators. 
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Appendix A: Some musings on PNLF and SBF 
distances 

One of the most intriguing issues in the extragalactic dis- 
tance scale is the ^ +0.3 mag a verage offset be tween the 
PNLF and SBF distance moduli (' Ciardullol[20T^ . Locally, 
there appears to be very little offset, but the discrepancy 
increases with distance, such that the PNLF method gives 
smaller mean distances for the Virgo and Fornax clus- 
ters, as well as a smal ler relative distance of Fornax with 
respect to Virgo (see IVilleeas et al.l 120101 ). Consistently, 
iFeldmeier et al.l ()20Q7D found that the value of Hq obtained 
from SNe la was 10% higher when the SN la distances were 
calibrated via PNLF distances, as compared to the Hq ob- 
tained by calibrating the SN la distances by either SBF or 
directly from Cepheid distances. 

As discussed in detail in Section fJH and ^ to reli- 
ably compare two or more distance indicators, it is of 
paramount importance to verify the consistency of the cal- 
ibrations (or calibrators) used. Both PNLF and SBF cali- 
brations are tied to the same primary indicator, the period- 
luminosity relation of C epheids, and to t he sam e dataset, 
i.e., the Cepheids from iFreedman et al.l ()2001t ). However, 
the zero points of present SBF calibrations are tied to the 
Cepheid distances obtaine d with metallicity-dependent PL 
relations (D pLZ hereafter: iMei et aLll2007t iBlakeslee et al.l 
120091 120101 this paper), while the standard PNLF calibra- 
tion relies on Cep heid distances with no dependence on 
meta l hcity (DpL '■ [ CiarduUo et al] 120021 : IFeldmeier et al.l 
I2OO7I: lCiarduTioll20m 

In th is appendix, we take the detailed discussion pre- 
sented in lCiardullo et al.l ()2002f ) - who found the ~ 0.3 mag 
mismatch between the two distance indicators - and ana- 
lyze the correction terms needed to homogenize the com- 
parison. 

First, as discussed above, iTonrv et al.l (|20Qll) distance 
moduli should be revised downward by 0.06 mag when us- 
ing the DpLz ■ For PNLF, instead, fciarduiy (\20T% re- 
ported a -1-0.07 mag correction to the PNLF distance mod- 
uli (0.07 mag brighter zeropoint, M*) when the Dplz dis- 
tances are used instead of I?pl . It is useful to note that 
the author also finds that the best-fit value to M* cali- 
brated against the RGB-Tip distances, i.e., independently 
from Cepheid distances, is again -1-0.07 mag brighter than 
the PNLF calibration obtained from DpL ■ 

The two corrections: a)— 0.06 mag for SBF distances, 
and b) -1-0.07 mag for PNLF, both deriving from the adop- 
tion of the DpLz , can justify 0.13 mag of the PNLF-SBF 
offset. 

Figure lA.ll shows the histogram of the differences be- 
tween PNLF and SBF distance moduli based on various 
SBF and PNLF distance and/or calibration. Panel a) in 
the figure shows the PNLF to SBF {m—M) difference. 
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Table 7. Quality control statistics for the exposures. 



Name 
- 




Background 
(counts) 


Background 
{mag / arcsec?) 


seeing 
(") 


Used Data 
- 






B band 








FORS, 


,1999-12-28T02:42:27.813.fits 


888.0 


22.25 


0.86 


yes 


FORS, 


,1999-12-28X02:54:22. 3'16.fits 


907.4 


22.22 


0.81 


yes 


FORS, 


,2000-01-09X02:41:46. 659.fits 


983.8 


22.13 


1.07 


yes 


FORS, 


,2000-01-09X02:52:43. 240.fits 


998.9 


22.12 


1.01 


yes 


FORS, 


,2000-01-12X01:11:59. 202.fits 


1556.3 


21.64 


1.06 


yes 


FORS, 


,2000-01-12X01:22:57.678.fits 


1527.3 


21.66 


1.12 


yes 


FORS, 


,2000-01-13X00:56:41. 525.fits 


2319.5 


21.20 


1.20 


yes 


FORS, 


,2000-01-13X01:07:40. 539.fits 


2236.8 


21.24 


1.14 


yes 


FORS, 


,2000-01-14X03:50:58. 405.fits 


1958.7 


21.39 


1.15 


yes 


FORS, 


,2000-01-14X04:01:58. 759.fits 


1580.1 


21.62 


1.19 


yes 


FORS, 


,2000-01-15X03:10:31. 420.fits 


7298.6 


19.96 


0.77 


no 


FORS, 


,2000-01-15X03:21:27.996.fits 


7231.2 


19.97 


0.86 


no 


FORS, 


,2000-01-17X03:28:21. 613.fits 


11877.6 


19.43 


1.36 


no 


FORS, 


,2000-01-17X03:39: 19. 186. fits 


12162.4 


19.40 


1.30 


no 


FORS, 


,2000-01-18X01:20:21. 447.fits 


12832.7 


19.35 


0.87 


no 


FORS, 


,2000-01-18X01:31:20. 684. fits 


12983.9 


19.33 


0.77 


no 


FORS, 


,2000-01-19X01:30:49. 601. fits 


20047.5 


18.86 


1.38 


no 


FORS, 


,2000-01-19X01:41:43. 505. fits 


21026.0 


18.81 


1.31 


no 


FORS, 


,2000-01-20X02:36:07. 201. fits 


32829.0 


18.33 


1.15 


no 


FORS, 


,2000-01-20X02:47:02. 044. fits 


37640.4 


18.18 


0.70 


no 


FORS, 


,2000-01-21X01:58: 13. 531. fits 


35260.4 


18.25 


0.70 


no 


FORS, 


,2000-01-21X02:09: 11. 551. fits 


35681.1 


18.24 


0.70 


no 






V band 








FORS, 


,2000-01-09X03:28:34.384.fits 


3261.5 


21.13 


0.91 


yes 


FORS, 


, 2000-0 1-09X03 : 39 : 32 . 352 .fits 


3333.3 


21.10 


1.10 


yes 


FORS, 


,2000-01-13X01:20:04. 617.fits 


10012.2 


19.91 


0.98 


yes 


FORS, 


,2000-01-13X01:31:04. 242.fits 


69.2 


19.93 


0.84 


yes 


FORS, 


,2000-01-13X01:39:32. 083.fits 


5141.9 


20.63 


1.03 


yes 


FORS, 


,2000-01-13X01:50:31. 241.fits 


5130.5 


20.64 


1.09 


yes 


FORS, 


,2000-01-17X04: 16: 19. 860.fits 


24234.6 


18.95 


1.51 


no 


FORS, 


,2000-01-17X04:27: 17.584.fits 


25234.8 


18.91 


1.59 


no 


FORS, 


,2000-01-18X02:06:59. 548.fits 


25642.0 


18.89 


0.76 


no 


FORS, 


,2000-01-18X02:17:56. 431.fits 


23947.3 


18.96 


0.93 


no 


FORS, 


,2000-01-19X01:07:22. 117.fits 


38819.6 


18.44 


1.21 


no 


FORS, 


,2000-01-19X01:18:19. 487.fits 


37909.1 


18.46 


1.28 


no 


FORS, 


,2000-01-20X01:52:46. 180.fits 


22531.1 


18.28 


1.25 


no 


FORS, 


,2000-01-20X01:58:44. 418.fits 


29233.8 


17.99 


1.36 


no 


FORS, 


,2000-01-20X02:05:45. 069.fits 


39325.9 


17.67 


0.95 


no 


FORS, 


,2000-01-20X02: 11:42. 821. fits 


38475.5 


17.70 


0.76 


no 


FORS, 


,2000-01-21X00:54:09. 577. fits 


6202.4 


17.93 


1.14 


no 


FORS, 


,2000-01-21X01:06:08. 681. fits 


5118.1 


18.14 


1.06 


no 


FORS, 


,2000-01-21X01:30:59. 576. fits 


16978.9 


18.14 


1.09 


no 


FORS, 


,2000-01-21X01:35: 17. 720. fits 


17047.7 


18.14 


1.11 


no 


FORS, 


,2000-01-21X01:39:35. 548. fits 


15737.2 


18.23 


0.98 


no 


FORS, 


,2000-01-21X01:43:50. 433. fits 


15366.4 


18.25 


1.03 


no 


FORS, 


,2000-01-21X01:48:07.248. fits 


15653.2 


18.23 


1.06 


no 


FORS, 


.2000-01-21T01:52:24.597.fits 


16583.5 


18.17 


1.08 


no 






/ band 








FORS, 


,1999-12-27X03:46:00. 987.fits 


9675.7 


19.07 


0.69 


yes 


FORS, 


,1999-12-27X03:57:55. 778.fits 


10469.2 


18.99 


0.78 


yes 


FORS, 


,2000-01-09X03:05:09. 924.fits 


8138.7 


19.26 


0.87 


yes 


FORS, 


,2000-01-09X03: 16:07.432.fits 


8755.6 


19.18 


0.82 


yes 


FORS, 


,2000-01-13X02:03: 14. 285.fits 


9061.0 


19.14 


0.93 


yes 


FORS, 


,2000-01-13X02:14:12. 265.fits 


9847.6 


19.05 


0.84 


yes 


FORS, 


,2000-01-15X03:34:22. 320.fits 


19878.0 


18.29 


0.79 


yes 


FORS, 


,2000-01-15X03:45:21. 020.fits 


20169.8 


18.27 


1.00 


yes 


FORS, 


,2000-01-17X03:52:07.180.fits 


18776.9 


18.35 


1.22 


no 


FORS, 


,2000-01-17X04:03:05. 404.fits 


19359.7 


18.32 


1.36 


no 


FORS, 


,2000-01-18X01:43:30. 890.fits 


19114.6 


18.33 


0.65 


no 


FORS, 


,2000-01-18X01:54:29. 531. fits 


19684.5 


18.30 


0.74 


no 


FORS, 


,2000-01-19X00:44:09. 869.fits 


27488.0 


17.91 


1.13 


no 


FORS, 


,2000-01-19X00:55:07.408.fits 


28624.7 


17.89 


1.09 


no 


FORS, 


,2000-01-21T02:22:08.776.fits 


16986.9 


17.71 


1.01 


no 


FORS, 


,2000-01-21T02:28:02.578.fits 


17131.9 


17.70 


1.05 


no 


FORS, 


,2000-01-21X02:33:59.741.fits 


16976.3 


17.71 


0.99 


no 


FORS, 


,2000-01-21X02:39:57.532.fits 


16960.7 


17.71 


1.08 


no 


FORS, 


,2000-02-04X02:06:34.213.fits 


9268.9 


19.12 


1.03 


yes 


FORS, 


,2000-02-04T02:17:28.971.fits 


8959.7 


19.16 


0.85 


yes 
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ApNLF-SBF, using; t he original sample of 28 galaxies by 
ICiardullo eraP (|2002[ ). with updated zero points for both 
distance indicators. The mean is Apnlf-sbf — —0.23 
mag, to be compared to Apnlp^sbf = —0.36 mag before 
zero-point correction. 

After the ICiardullo etaD (|2002l) paper, few PNLF dis- 
tances have been obtained for galaxies with SBF measure- 
ments. Panel b) in Figure lA.ll shows the PNLF to SBF 
difference for a total of 33 galaxies, including the dis- 
tances obtained after 2002. The comparison shown in panel 
c) is obtained using the recent SBF distances from the 
ACSV CS and ACSFCS s urveys (when available) in place 
of the iTonrv et al.l ()200lD distances. Finally, panel d) uses 
the same SBF and PNL F distances of panel c) except that 
for the galaxies with old iTonrv et al.l (|200l[ ) corrected dis- 
tances w e include the further " Q-correction" term using eq. 
Al from lBlakeslee et al.l (|2010f ). Table lA. II presents the av- 
erage and median differences for all assumptions shown in 
Figure EH 

It is wo rth mentioning that, according to 

ICiardullo et"all (|2002[) . to properly compare PNLF 
and SBF, the latter distance moduli should be further 
reduced by 0.04 mag. If one includes this correction term, 
all differences reported in Table lA.ll becomes smaller, 
with the best sample (i.e., d in the table) providing 
ApNLF-SBF — —0.14 ± 0.06 mag, and a median of —0.22 
mag. 

Taking into account all the corrections described above, 
the offset between SBF and PNLF distances is reduced to 
half of the original estimate, once the proper calibrations 
and recent /updated distances are used for both indicators. 
However, even in the best case there is a non-negligible 
^0.2 mag offset between SBF and PNLF that implies 
~ 10% larger SBF distances. Again, this difference occurs 
mainly beyond ^ 10 Mpc, and is similar to the 10% larger 
Hq obtained when calibrating SNe la via PNLF instead of 
Cepheids (iFeldmeier et al. 200 7j) . 

Even though the data presented in Table lA.ll (with the 
possible further -1-0.04 mag improvement cited above) sug- 
gest that the best average difference is statistically consis- 
tent with zero - especially taking into account the system- 
atic uncertainties, not considered in this comparison - the 
scatter wc find is larger, or nearly equal to the squared sum 
of the estimated internal scatters of both indicators. Taken 
at face value, this result either means that the internal scat- 
ter of one or both indicators is underestimated, or that a 
real systematic offset exists between the two. 

As another test, to further check the latter concern, 
we analyzed the PNLF to SBF offset by considering late- 
type and early-type galaxies separately. Using sample d in 
Table rSm we find the following differences: A^°^^j^p_ggp = 
—0.23 ±0.05 mag (median —0.30 mag) based on the data of 
24 galaxies, and Ajif^'ip^ggp — 0.04 ± 0.17 mag (median 
—0.02 mag) for the remaining 9 galaxies. The result shows 
that for the class of galaxies used to derive SBF and PNLF 
zero points, i.e., the late types hosting Cepheids, there is 
no statistically significant offset between the two indicators. 
On the other hand, for early-type galaxies the offset is large, 
and statistically inconsistent with zero. Related to this, one 
fundamental difference between the two indicators is that 
absolute SBF magnitudes are "corrected" for galaxy stel- 
lar content, i.e., the difference between early- and late- type 
galaxies is taken into account with SBF, while the PNLF 



distances are based on the constancy of M* for both types 
of galaxy. 

The calibration of SBF magnitudes, and its depen- 
dence on galaxy type, has been analyzed in detail over 
a grou nd-based sample of ^^300 galaxies by iTonrv et al.l 
(1200 Ih. and more r e cently from HST data of ~ 150 galaxies 
bv iBlakeslee et al.l (|2009[ l2010D . The derived SBF calibra- 
tions, as is well-known, include a colour- dependent term , 
which, as also shown by SSP models (IWorthevI Il993bl : 
ICantiello et al.l [20031 : iRaimondo et al.ll2005f) . is basicallv a 
metallicity correction term. In optical bands, this correc- 
tion term implies fainter SBF magnitudes for redder/more 
metal-rich systems. 

In contrast, the M* calibration to PNLF does not in- 
clude any metallicity dependent term for bright galaxies. A 
de pendence of the PNLF M* to metallici ty has been found 
by ICiardullo etlD (|2002l) and ICiardullol ((20JJ), however 
the authors conclude that such dependence is relevant only 
in small, metal-poor systems. 

In specting the open circles in Figure 5 of ICiardullol 
(|2012l open circles mark the data obtained from DpLz dis- 
tances) , reported in Figure IA.21 one can see that a) the 
trend of M* with metallicity has the opposite sign with 
respect to SBF, meaning that M* gets brighter for more 
metal-rich systems, and that b) some residual correlation 
of M* with metallicity also appears in the high metallicity 
regime. 

This strongly suggests that there may be some unac- 
counted for residual dependence of the PNLF calibration on 
the metallicity, presently unquantified because of the rela- 
tively limited sample - though detectable even in present 
datasets (Fig. lA.2| ). If so, one possible explanation for the ir- 
reducible ~ —0.2 mag offset between PNLF and SBF is that 
it may be due to a bias introduced by data from early-type 
galaxies, which are intrinsically on average more massive 
and more metal-rich than late types and, therefore, would 
typically have larger distances with respect to the current 
metallicity-independent calibration of the PNLF method. 



Appendix B: Some notes on Type la SNe distances 
from Stritzinger et al. (2010) 

To understand the possible causes of the difference between 
our and StrlO distance we must recall that, as for SBF, 
in order to calibrate SNe la light curves one must rely on 
sources at known distance and/or with well-known intrinsic 
properties, and then standardize the absolute magnitude of 
the SN la (e.g., [Phillips 1993). 

StrlO derives the distance to NGC 1316 using three dif- 
ferent methods: the EBV, the Tripp method, and the near- 
IR light-curves. The authors also use the MLCS2k2 method 
on SN 2006mr, obtaining a distance that is 50% further than 
the average they derived from the normal events. 

In the following we discuss each one of the three meth- 
ods used by StrlO, trying to highlight the possible causes 
leading to the observed difference. 

B.l. The Colour Excess Method 



StrlO adopted the calibration from iBurns et al.l (j201ll AJ 
in press at the time StrlO was published) to obtain the dis- 
tance moduli with the "EBV" model in their fitting package 
SNooPy. For this model they adopt the recommended cal- 
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Table A.l. SBF to PNLF comparison. 



Sample 


Number of galaxies 


ApNLF-SBF {rms) 
(mag) 


Median difference 
(mag) 


Original'^' 


28 


-0.36 (0.31) 


-0.36 


al^l 


28 


-0.23 (0.31) 


-0.23 


b[3] 


33 


-0.21 (0.29) 


-0.23 


cW 


33 


-0.22 (0.30) 


-0.29 


d[5l 


33 


-0.18 (0.34) 


-0.26 



Ciardullo et al. 



Ciardullo et al 



(2002) with SBF distances from lTonrv et al.l l|2001l '). 

( 20021 ) PNLF with revised zero points based on metallicity dependent PL relations 



] Original PNLF sa mple by 
j iTonrvet all (|2001h SBF & 

f or Cepheids ( s ee tex t). 

] iTonrv et all (2001') SBF & ICiarduUo eF aT (2002* ) PNLF plus more recent PNLF dis tances fr om Fcldm eier et al.1 l|200' 
NGC 1380, NGC 4526 and r evised NGC 1316 ), Sani bhus et al.l (2006, NGC 4697), Tcodoresc^ueFaD ([2005, 2010' NGC 1344 and 
NGC 821, respectively) and lHerrmann et al.l (|2008l . NGC 4376). For both distance indicators revised zero points jaxe used, as in 
sample b except for NGC 4697, NGC 1344 and NGC 821 who se PNLF dist ance is independent from iCiarduUo et al.l calibration. 
] Updated SBF distances from ACSVCS (|C6te et al.ll2004l ) & ACSFCS (| Jordan et al.l 120071) (when available). Old SBF and 

PNLF distances as in sample b. 

] As sample c except that SBF with Q-corrected SBF distances (see text) are used for the old I Tour v et al.l (|200ll ) distance 
moduli. 
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Fig. A.l. Histograms of the difference between PNLF 
and SBF distance moduli. Each panel shows different 
comparison/calibration samples. Panel a): Original sam- 
ple of 28 ga laxies with PNLF and SBF di stances from 
ICiardullo et al. (2002) and lTonrv eraD (l200lt) respectively, 
with revised Cepheids calibration. Dotted/dashed vertical 
lines show the mean/median of the difference. Panel b) : 
As in a) but 5 more galaxies w ith PNLF measurements 
made after ICiardullo et al.l (|2002f ) are added to the sample. 
Panel c): As h), but more recent ACSFCS and ACSVCS 
SBF distances are used when av ailable. Panel d ) -. As c), 
but the cor rection from e g. Al bv lBlakeslee et al.1 ()2010[ ) is 
included to iTonrv et al.l (|200lh distances. 



ibration from iBurns et al.l (l201lh. which us e s a su b-sample 
of unreddened SNe la from iFolatelh et all (|2010D. exclud- 
ing fast declining objects, as SN2006mr. Folatelli et al.l 
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Fig. A. 2. Figure 5 from ICiardullol (j2012|) showing the sole 
data with metallicity corrected Cepheids distances. [See 
electronic version of the Journal for a colour version of 
the figure.] 



()2010[ ). in turn, calibrate their dataset using 26 SNe la at 
z > 0.01, whose distances are based on Bubble's law assum- 
ing Ho = 72 km s"^ Mpc"\ and three SNe la at z < 0.01 
with direct distance measurements. The three nearest ob- 
jects include SN2006mr, whose host's distance is assumed 
(m — M) = 31.59 =b 0.0 8 mag from the SBF measurements 
bv lCantiello et al.l (pOOTa ). 

The SNooPy/EBV method relies (also) on an estimate 
of the internal reddening around the SN la. StrlO find that 
all four SNe la in NGC 1316 have negligible internal ex- 
tinction. The authors warn about some complications in 
the interpretation of the data. With respect to spectro- 
scopic analysis of Na I D absorption in the spectra of 
SN2006dd and SN2006mr, these authors state that "the 
very strong Na I D absorption observed in SNe la 2006dd 
and SN 2006mr is totally inconsistent with the low host 
galaxy reddening we derive from the light curve observa- 
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tions" (StrlO, Section 4), as the colour evolution of the 
two SNe la closely resembles that of unreddened SNe la. 
However, while there is a general observational agreement 
on higher colour excess corresponding to higher Na I D 
equivalent width (EW), this correlation is tight in high- 
resolution spectra, b ut the scatter increases substantially 
at lower re solution ("Blo ndin et al.l 120091 : iPoznanski et alJ 
1201 ll [20T2l ). with considerable confusion due to the blend- 
ing of the Na I D do ublet. For instance, using Figure 5 in 
iBlondin etHI (|2009D . at EW~ 1.5 (similar to that obtained 
by StrlO for NGC 1316) E{B - V)host ranges from - 0.1 
to 1.5 mag. As a further complication, besides the possi- 
ble line blending, the line profiles for both SNe show clear 
evidence for structure, interpreted by StrlO as evidence 
of the presence of two unresolved sodium components. 
Nevertheless, coupling the EW measurements from StrlO 
(their Table 10, obtained from the author s' highest disper- 
sion s pectra), with equations (7)-(8) from IPoznanski et al.l 
(|2012D . we obtain E{B - V)host > 0.45 and > 0.15 mag for 
SN 2006dd and SN 2006mr, respectively (with mean values 
~ 1.4 and ~ 0.6 mag). 

Clearly, as also highlighted by StrlO, the strong Na I 
D absorption associated with no internal reddening might 
also be indicative of a non-standard gas-to-dust ratio for 
Fornax A, which given the merger history of this galaxy 
does not seem unreasonable. 

A 6" X 6" zoom of the regions around the four SNe la 
is shown in Figure IB. II In the upper part of the fig- 
ure we show the HST UVIS/F336VF, ACS/F435VK and 
ACS/F555W residual frames of the three SNe la lo- 
cated within the frames analyzed in this work (SN 1981D, 
SN2006dd, SN2006mr). For sake of completeness, we have 
obtained archival g, r, and z-band Gemini/GMOS-S data of 
the region around the SN 1980N, shown in the lower panels 
of Figure lB.ll As evidenced in the figure, there are undeni- 
able patterns of dust near the position of the 2006 SNe la, 
though one cannot decide whether the SNe are behind, in 
front of, or within such dust lanes. 

The comparison with previous estimates of internal ex- 
tinction fr om the lit erature show that for SN 1980N and 
SN 1981D. lJha et all (12007 ) find E{B-V)host values larger 
than StrlO (three times larger in the case of SN 1981D) 
but in agreement within uncertainties with the SNooPy 
fits. The agreement gets worse if the extinctions derived 
by StrlO from near-IR data are taken into account. Note, 
however, that StrlO corrected the optical and near-IR pho- 
tometry for host galaxy contamination. Such a correction, 
although negligibl e for the ca se of SN 1980N, does not seem 
to be discussed bv lJha et al.1 

For both the two mos t rece nt SNe la, 2006dd and 
2006mr,|M aoz fc Mannuccil ()2008[ ) estimate an internal ex- 
tinction of -0.08 mag. In the case of SN 2006dd, SNooPy 
provides E{B—V)host = 0.043 ± 0.008 mag; while, as men- 
tioned above, SNooPy cannot be used for fitting the light- 
curves of the fast declining SN 2006mr. 

In conclusion, the three SNe la used by StrlO to get the 
best estimate of (m— M) could be controversial in terms of 
internal extinction, affecting both the estimate of (m— M) 
with the EBV method, and the associated uncertainties. 

B.2. The Tripp Method 

The second metho d adopted by StrlO is based on the two- 
parameter model of lTripd (jl998[ ) which, differently from the 



EBV method, can also be applied to fast de clining SNe la. 
The calibrating sample is again taken from iFolatelli et ah! 
(2010), and SN2006mr is omitted in the re-computed cali- 
bration relations, to avoid circularity. 

The distance to the three normal SNe la with this 
method is consistent with the estimates based on the 
SNooPy/EBV method. In contrast, the (m—M) obtained 
with the data of SN2006mr is — -1-0.5 mag larger than the 
average of the other three SNe la. 

One interesting point to note is that if one uses a dis- 
tance modulus of ~31.2 mag to NGC 1316, as derived by 
StrlO from the three normal SNe l a, then the dat a -point of 
SN2006mr placed in Figure 16 of IFolatelli et al.l ()2010t) is 
more than -1-0.5 mag off the linear relation drawn by the au- 
thors, with a scatter to the relation much larger than the 
rms reported in the cited figure. Certainly, changing the 
dist ance modulus of one of the calibrating data-points in 
the IFolatelli et al.l sample implies changing the linear cali- 
bration relation shown in the cited figure, and possibly re- 
duces the offset between data and fit. In any case, though, 
using (to — M) ^ 31.2 entails considerably increasing the 
scatter in the calibrating sample of the Canergie Supernova 
Project, with NGC 1316 being one of the nearest objects 
and also the host galaxy with the largest scatter. 

B.3. The near-IR method 

The last method used by StrlO is base d on near-IR light- 
curves of SNe la, calibrated using Kris ciunas et al.l ()2009() 
abso lute near-IR peak magnitudes without NGC 1316 data. 
The E risciunas et al.l calibration of near-IR peak mag- 
nitudes adopts new observations of S Ne la, and data 
previously published by the same team ([Krisciunas et aTj 
2004a, b). For the nearby galaxies, the authors adopted dis- 
tances based on either SBF or Cepheids . From a careful 
reading of the cited papers, we find that iKrisciunas et al.l 
()2004aD obtained the JHK calibration from 16 SNe la. 
For three nearby galaxies, NGC 1316, NGC 4526, and 
NGC 5128, t he au thors adopt the SBF distance from 
lAihar et"al] (|2001l based on Dplz distances), while for 
NGC 4536 and NGC 3368 the L»pl distances from Cepheids 
is use d. Both SBF and Cepe hids distances are based on the 
same 'Freedman et al." (2001') calibrating sample. 

fKrisciun as et a l. (2004b,) extended the sample of SNe la 
with well-sampled near-IR light-curves to about 20 objects. 
The authors added two more supernova-host galaxies with 
SBF distances - NGC 4374 and NGC 3190 - to the pre- 
vious list of nea r by gal axies. However, in c ontrast with 
IKrisciunas et al.l (|2004al ). they adopted the iTonrv et al.l 
()200lD distance moduli, which are based on Dpl ; that 
is, they are 0.06 m ag larger than the (m—M) reported by 
lAihar et al.l (|200lh . 

Finally, the most rec ent calibrat i on of near-IR light- 
curves of SNe la bv .Krisciunas et all (|2009D . used the dis- 
tance modul i for ri earby galaxies fron i the SBF survey by 
iTonrv et all (|200lh with the revised Ijensen et al." (2003() 
zero points (+0.16 mag with respect to Tonry ot al. 2001). 
More specifically, in the new list of ~25 supernova-host 
galaxies, the authors added three new objects, NGC 936, 
NGC 1 201 and NGC 137 1, with SBF-based distances. In ad- 
dition, 'Krisciunas et al.' adopted a revised SBF distance to 
NGC 5128 from Jensen ct al. (2003). However, the authors 
corrected the (to— M) by +0.16 m ag, which is the difference 
between the iTonrv et all (|200lh and iJensen et af] (|2003D 
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SN 1980 N 

Fig. B.l. Zoom in 6 x Qarcsec^ of the residual frames centered on the positions of the four SNe la, as labeled. Black 
circles mark the position of the respective SNa. UVIS/i^336W^, ACS/i^435W^ and ACS/F555W^ frames (left, middle and 
right panels, respectively) are shown for the SN1981D, SN2006dd and SN2006dd. The panels showing SN 1980N are g, 
r, and i observations taken from the Gemini Telescope archive. 



calibratio ns, and not by +0.1 m ag, i.e., the differenc e be- 
tween the lAihar eFaLl ()200lD and lJensen et al.l (|2003D cal- 
ibrations. 

Hence, the calibration of the absolute magnitudes in 
near-IR bands used by StrlO relies on a sample of '--^25 
SNe la, with nine nearby calibrators (after excluding 
NGC 1316) having distances based on non-homogeneous 
calibrations . In Table IB . 1 1 we summarize the distance mod- 
uli used by iKrisciunas et al.l to calibrate the maxima of 
SNe la light-curves, and the ones revised by us in order to 
be a) internally homogeneous, and 6) consistent with the 
Cepheid distances used in this work. As shown in the table, 
the revised distance moduli are on average 0.06 mag larger 
than the ones used for the original SNe la near-IR cali- 
bration. If one simply takes the average of these numbers, 
adding a sample of 15 objects more (25 total SNe la minus 
the nine nearby objects and NGC 1316) where no shift has 
to be applied, the correction to the absolute magnitudes in 
Table 9 of StrlO is ~ —0.02 ma g (a correction that sho uld 
be applied to the calibration bv IKrisciunas et 311120091) . 

Although a 4-0.02 mag shift in StrlO distance moduli 
goes in the direction of reducing the mismatch between our 
and StrlO distances, the amplitude of the correction is negli- 



gible. However, it suggests that the uncertainties associated 
with the near-IR calibrations might be underestimated. 

In conclusion, the analysis presented in this Appendix high- 
lights two main issues: i) the homogeneity of the calibra- 
tors used, and ii) the estimate of internal extinctions for 
SNe la (although one should not forget that the quality 
of the SN 1980N and SN 1981D data is lower than others). 
While the first issue listed works in the direction of reduc- 
ing the difference between our and StrlO distance moduli, 
both issues imply an increase of the present levels of statis- 
tical/systematic uncertainties on SNe la distances. 
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Table B.l. Nearby galaxy sample used for near-IR calibration: original and revised distances. 



Galaxy 


{m — M)near-IR 


Reference 


Method 


{m — M) Revised 


A(Revised — Orig.) 


NGC 4526 


31.08 


Krisciunas^t_al^ (2004a) 


SBF 


31.08 


+0.00 


NGC 4536 


30.80 


Krisciunas et al. f2004a) 


Cepheids 


30.87 


+0.07 


NGC 3368 


29.97 


Krisciunas et al. f2004a") 


Cepheids 


30.11 


+0.14 


NGC 4374 


31.32 


Krisciunas et al. f2004b) 


SBF 


31.26 


-0.06 


NGC 3190(NGC 3226) 


31.86 


Krisciunas et al. f2004b) 


SBF 


31.80 


-0.06 


NGC 936 


31.65 


Krisciunas et al. (2009) 


SBF 


31.75 


+0.10 


NGC 1201 


31.37 


JCrisciunas et_al^ (2009) 


SBF 


31.47 


+0.10 


NGC 1371(Eridanus group) 


31.84 


Krisciunas et aL (2009) 


SBF 


31.94 


+0.10 


NGC 5128 


27.90 


Krisciunas et aL (2009) 


SBF 


28.06 


+0.16 



Average (median) correction on the sample of nearby galaxies 0.06+0.03 (0.085) 

Average correction on the complete sample of galaxies (assume zero for distant galaxies) 0.02+0.01 
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